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ABSTRACT 

Context. Accurate mass, radius, and abundance determinations from binaries provide important information on stellar 
evolution, fundamental to central fields in modern astrophysics and cosmology. 

Aims. We aim to determine absolute dimensions and abundances for the three F-type main-sequence detached eclipsing 
binaries AD Boo, VZ Hya, and WZ Oph and to perform a detailed comparison with results from recent stellar evolutio- 
nary models. 

Methods, uvby light curves and uvby[3 standard photometry were obtained with the Stromgren Automatic Telescope 
at ESO, La Silla, radial velocity observations at CfA facilities, and supplementary high-resolution spectra with ESO's 
FEROS spectrograph. State-of-the-art methods were applied for the analyses: the EBOP and Wilson- Devinney binary 
models, two-dimensional cross-correlation and disentangling, and the VWA abundance analysis tool. 
Results. Masses and radii that are precise to 0.5-0.7% and 0.4-0.9%, respectively, have been established for the com- 
ponents, which span the ranges of 1.1 to 1.4 Mq and 1.1 to 1.6 Rq. The [Fe/H] abundances are from —0.27 to +0.10, 
with uncertainties between 0.07 and 0.15 dex. We find indications of a slight a-element overabundance of [a/Fe]~ +0.1 
for WZ Oph. The secondary component of AD Boo and both components of WZ Oph appear to be slightly active. 
Yale-Yonsai and Victoria-Regina evolutionary models fit the components of AD Boo and VZ Hya almost equally well, 
assuming coeval formation, at ages of about 1.75/1.50 Gyr (AD Boo) and 1.25/1.00 Gyr (VZHya). BaSTI models, 
however, predict somewhat different ages for the primary and secondary components. For WZ Oph, the models from all 
three grids are significantly hotter than observed. A low He content, decreased envelope convection coupled with surface 
activity, and/or higher interstellar absorption would remove the discrepancy, but its cause has not been definitively 
identified. 

Conclusions. We have demonstrated the power of testing and comparing recent stellar evolutionary models using eclipsing 
binaries, provided their abundances are known. The strongest limitations and challenges are set by T a a and interstellar 
absorption determinations, and by their effects on and correlation with abundance results. 

Key words. Stars: evolution - Stars: fundamental parameters - Stars: binaries: close - Stars: binaries: eclipsing - 
Techniques: photometry - Techniques: spectroscopic 



1. Introduction 

Accurate stellar mass and radius data from eclipsing bi- 
nary systems are valuable empirical test data for stellar 
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evolution models, basically because they are free of any 
scale dependent calibrations. Additional tests become pos- 
sible when the effective temperatures for the components 
can also be determined, most often from well-calibrated 
photometry. The binaries then also serve as primary dis- 
tance indicators, e.g. for nearby stars, stellar clusters, and 
Local Group galaxies, and they define the empirical mass- 
luminosity relations over a broad mass range. However, the 
most stringent model tests require that spectroscopic ele- 
ment abundances are available as well, minimising the num- 
ber of free parameters in the comparison with theory (e.g. 
Andersen ll99ip . Unfortunately, the number of binary com- 
ponents with complete and accurate data is still small. 

In our efforts to provide such complete data for a larger 
sample of eclipsing binary systems, F-G type main-sequence 
stars have been given high priority for two reasons. First, 
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their spectra offer favourable conditions for accurate radial 
velocity determinations and abundance analyses. Second, 
such stars cover the mass range within which convective 
cores begin to develop and affect the observational deter- 
mination of stellar ages, in particular due to difficulties 
and shortcomings related to the physics of core overshoot. 
Moreover, these ages cover the range of interest in studies 
of the evolution of the Galactic disk (e.g. Holmberg et al. 
120071) . 

In this paper we present new analyses of the three 
detached F-type double-lined eclipsing binaries AD Boo, 
VZ Hya, and WZ Oph. From earlier studies they are known 
to have masses in the range 1.15-1.41 Mq and radii cor- 
responding to the lower half of the main-sequence band, 
but better accuracy and - especially - abundance data are 
needed for meaningful tests of state-of-the art stellar mo- 
dels. We have therefore obtained new radial velocity ob- 
servations, light curves, standard photometry, and high- 
resolution spectra covering a wide wavelength range, en- 
abling us to determine all relevant parameters of the sys- 
tems with high precision. 

As pointed out by Popper (|1998ap , and also seen in 
Fig. [Dlj, AD Boo is of particular interest for having com- 
ponents with the greatest differences in mass and radius 
among well-studied binaries in its mass range. Hence, it 
may provide stringent tests of stellar models, but the 
masses obtained in two recent analyses of the system dis- 
agree by slightly more than their stated errors (Lacy[1997a; 
Popper 1998a]). Moreover, in both studies the radii were 
based on reanalyses of the relatively old B, V light curves 
(normal points) by Zhai et al. (|1982p . which have an accu- 
racy of only about 0.02 mag and show out-of-eclipse varia- 
tions by about 0.075 mag. 

Previous determinations of the dimensions of VZ Hya 
and WZ Oph were also based on old material fPopper ["1965p 
and therefore had an accuracy of only 2-4%; see Popper 
(j!980|) . As in the case of AD Boo, the components of VZ Hya 
are fairly different, whereas those of WZOph are nearly 
identical. 

The paper is structured as follows: First, Sect. [2] de- 
scribes our observations and analysis techniques adopted. 
The spectroscopic and photometric analyses of the indivi- 
dual systems are then discussed in detail in Sect. 03 The re- 
sulting new, accurate absolute dimensions and abundances 
are presented in Sect. [4] and compared with three recent 
grids of stellar evolutionary models in Sect. [5] 

Throughout the paper, the component eclipsed at the 
deeper eclipse at phase 0.0 is referred to as the primary (p) 
component, the other stars as the secondary (s). 



2. Observations and analysis techniques 

2.1. Spectroscopic observations 

The radial velocity observations for the three binaries 
were conducted at the Harvard- Smithsonian Center for 
Astrophysics (CfA) with nearly identical echelle spec- 
trographs on the 1.5-m Wyeth reflector at the Oak 
Ridge Observatory (Harvard, Massachusetts), the 1.5-m 
Tillinghast reflector at the F. L. Whipple Observatory (Mt. 
Hopkins, Arizona), and the Multiple Mirror Telescope (also 
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Fig. 1. Main-sequence eclipsing binaries with both com- 
ponents in the 1.1-1.5 Mq range having masses and radii 
accurate to 2% or better. The nine systems shown are, in or- 
der of increasing mass of the primary component: EWOri, 
UXMen, HSHya, V505Per, ITCas, V1143Cyg, CD Tau, 
DMVir, and BWAqr. WZOph, VZHya, and AD Boo, are 
shown with star symbols. The full drawn line is the Y 2 0.3 
Gyr isochrone for Z = 0.0181 (Demarque et al. I2004|) . 



on Mt. Hopkins) prior to its conversion to a monolithic mir- 
ror. A single echelle order 45 A wide was recorded with in- 
tensified photon-counting Reticon detectors, at a central 
wavelength of approximately 5187 A, which includes the 
Mg I b triplet. The resolving power of these instruments 
is A/AA w 35 000, and the signal-to-noise (S/N) ratios 
achieved per resolution element of 8.5 kms -1 range from 
10 to 43 for AD Boo, 15 to 39 for VZ Hya, and 10 to 25 for 
WZ Oph. As seen in Figures [H [9l and [TU the observations 
cover the orbits well. 

For detailed abundance studies, we have obtained high- 
resolution spectra at higher S/N ratios with the FEROS 
fiber echelle spectrograph at the ESO 1.52-m and 2.2-m 
telescopes at La Silla, Chile (Kaufer et al. [19991 [2000]) . The 
spectrograph, which resides in a temperature-controlled 
room, covers without interruption the spectral region from 
the Balmer jump to 8700 A, at a constant velocity resolu- 
tion of 2.7 kms -1 per pixel (A/AA = 48 000). An observing 
log is given in Table [1] 

A modified version of the MIDAS FEROS package, pre- 
pared by H. Hensberge, was used for the basic reduction of 
the FEROS spectraQj. Compared to the standard package, 
background removal, definition and extraction of the indi- 
vidual orders, and wavelength calibration are significantly 
improved. The observations were reduced night by night 
using calibration exposures (ThAr and flat field) obtained 
during the afternoon/morning. Standard (rather than op- 
timal) extraction was applied, and no order merging was 
attempted. Typically, a standard error of the wavelength 
calibration of 0.002-0.003 A was obtained. Subsequently, 
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Table 1. Log of the FEROS observations of AD Boo, 
VZ Hya, and WZ Oph. HJD refers to mid-exposure. The 
exposure time t cxp is given in seconds. The signal- 
to-noise ratio per resolution element (S/N) of individ- 
ual exposures was measured around 6070 A. Observers: 
H = Heidelberg/ Copenhagen guaranteed time; JVC = 
J.V. Clausen; JS = J. Southworth. 
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dedicated IDL0 programs were applied to remove cosmic 
ray events and other defects, and for normalization of the 
individual orders. For each order, only the central part with 
acceptable signal-to-noise ratios was kept for further ana- 
lysis. A short region of sample FEROS spectra is shown in 
Fig.! 

2.2. Radial velocity determinations 

Radial velocities were obtained from the CfA spectra using 
the two-dimensional cross-correlation algorithm TODCOR 
(Zucker & Mazeh [1994p . Templates were selected from a 
library of synthetic spectra based on model atmospheres 
by R. L. Kurucz (see also Nordstrom et al. 1994; Latham 
et al. I2002p . calculated over a wide range in effective tem- 
perature (T c ff), rotational velocity (vsini), surface grav- 
ity (log g), and metallicity ([M/H]). The templates produc- 
ing the best match to each component of the binary were 
identified by running extensive grids of cross-correlations 
and selecting the combination giving the highest correla- 
tion value averaged over all exposures, weighted according 
to the strength of each spectrum (e.g. Torres et al. [2002). 
These grids were in general run in temperature and rota- 
tional velocity only, since those are the parameters that 
affect the velocities the most. Surface gravities were fixed 
according to preliminary estimates. If needed, grids for dif- 
ferent assumed heavy element abundances were compared, 
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leading indirectly to a rough estimate of [M/H] for the sy- 
stem, which was in turn compared with the results from 
the detailed abundance analysis. 

The stability of the zero-point of our velocity system was 
monitored by means of exposures of the dusk and dawn sky, 
and small run-to-run corrections were applied as described 
by Latham 0992). In addition to the radial velocities we 
derived the light ratio at the mean wavelength of our ob- 
servations (5187 A), following Zucker & Mazeh (fl994j) . 

Because of the narrow wavelength coverage of the CfA 
spectra, the radial velocities were checked for systematic 
errors that might occur as a result of lines of the stars mo- 
ving in and out of the spectral window with orbital phase 
(Latham et al. I1996p . Numerical simulations were carried 
out as described in more detail by Torres et al. (|1997p . and 
corrections based on them were applied to the raw veloci- 
ties. These corrections are included in the measurements 
reported in Tables |A~T1 IA~2l andlA~3l 

2.3. Photometric observations 

uvby light curves for AD Boo, VZ Hya, and WZ Oph were 
obtained with the Stromgren Automatic Telescope (SAT) 
at ESO, La Silla during several campaigns between March 
1988 and March 1997. They contain 652, 1180, and 697 
points in each band, respectively, and cover all phases well. 
The accuracy per point is 4-5 mmag in vby and 6-7 mmag 
in u, far better than previous photometry of these systems. 
The light curves are kept in the stable uvby instrumental 
system of the SAT in order not to introduce additional 
scatter from transformation to the standard uvby system. 

Improved linear ephemerides have been calculated for 
all three systems from the new and numerous published 
times of minima. No signs of period changes were detected, 
and the orbits were assumed to be circular since there is no 
evidence for eccentricity, either in previous data or in our 
own. 

In addition to the light curves, homogeneous standard 
uvby [3 indices have been established from SAT observations 
on dedicated nights in which the binaries (and comparison 
stars) were observed together with a large sample of uvby 
and P standard stars. 

We refer to Clausen et al. PUUSI hereafter CVG08) for 
further details on the photometry and the new ephemerides. 

2.4. Photometric analyses 

The choice of binary model for the light curve analysis 
depends strongly on the relative radii of the components 
and the irradiation between them. In our case, the sy- 
stems are well-detached with relative radii between 0.10 
and 0.17, the components are little deformed (oblateness 
between 0.001 and 0.007), and reflection effects are small. 
Therefore, the simple Nelson-Davis-Etzel model (Nelson & 
Davis [1972 Etzel fTMTl Martvnov [T973)) . which represents 
the deformed stars as biaxial ellipsoids and applies a sim- 
ple bolometric reflection model, is expected to be perfectly 
adequate. We have adopted the corresponding EBOP code 
(Popper & Etze l | 1981"|) . supplemented by an extended ver- 
sion JKTEBOFQ including Monte Carlo simulations (see 
e.g. Southworth et al. I2004al 12004b)) . for the analyses of 
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Fig. 2. FEROS spectra (black line) of AD Boo (phase 0.17), VZ Hya (phase 0.28), and WZ Oph (phase 0.33) in the region 
of the Fel lines at 6056.0 and 6065.5 A, shifted to laboratory wavelengths for the primary components. Synthetic spectra 
have been included: primary components (dashed blue); secondary components (dashed green); combined (continuous 
red line). 



the three systems and the assignment of realistic errors for 
the photometric elements. 

Since the accuracy of the light curves is high, we have 
also checked for model-dependent systematic errors in the 
photometric elements by performing parallel analyses with 
two more advanced binary models, WINK (Wood 1972) and 
the widely used Wilson-Devinney (WD) model (Wilson & 
Devinney [Uni Wilson ITM31 H554]) . WINK uses tri-axial 
ellipsoids, whereas WD is based on Roche geometry. Our 
WINK version is modified and extended as described by 
Vaz <TT984l 11986b . Vaz k Nordlund (fl985l) . and Nordlund 
& Vaz (1990). Our WD code is a revised and extended 
version of the 1994 code by Wilson; see Vaz et al. (|1995p and 
Casey et al. (jl998). In both codes we use the predictions of 
ATLAS9 model atmospheres (Kurucz ll993p to connect the 
flux ratios of the two components in different passbands. 

In their publicly available versions, EBOP, WINK, and 
WD carry out a least-squares optimization of the pa- 
rameters by means of the classical differential corrections 
method. In JKTEBOP, this has been replaced by the 
Levenberg-Marquardt minimization algorithm (MRQMIN: 
Press et al. HMD). 

Throughout this paper we have analysed the uvby light 
curves independently, and for each band identical weights 
have been assigned to all observations. A photometric scale 
factor (the magnitude at quadrature) was always included 
as an adjustable parameter, and the phase of primary 
eclipse was allowed to shift from 0.0. The mass ratio be- 
tween the components was kept at the spectroscopic value, 
and synchronous rotation was assumed. Gravity darkening 
coefficients/exponents corresponding to convective atmo- 



spheres were applied. The bolometric reflection albedo co- 
efficients were fixed at 0.5 in the WINK and WD analyses, 
again due to convection; for EBOP the simple bolometric 
reflection model was used. 

Linear limb-darkening coefficients were either assigned 
from theoretical calculations (van Hamme 119931 Diaz- 
Cordoves et al. HME[ Claret l2"UUDj) according to the ef- 
fective temperatures, surface gravities, and abundances, or 
included as adjustable parameters. For all components of 
the three systems, the linear coefficients by Claret (2000) 
are about 0.10 higher than those by van Hamme |1993) 
and about 0.03 higher than those by Diaz-Cordoves et 
al. (pL995). In general, a change in the coefficients by 0.10 
causes systematic differences between the optimum theore- 
tical light curves of just a few mmag through the eclipses, 
and the overall quality of the fits to the observations is 
nearly identical. If adjusted, the coefficients for the com- 
ponents of VZ Hya and WZ Oph turn out closest to the 
values by van Hamme; for AD Boo the coefficients either 
become unrealistic or convergence is not achieved at all. 
In addition, a few WD analyses with the non-linear law 
by Klinglesmith & Sobieski (|1970p and coefficients by van 
Hamme (1993J were included. No improvements in the light 
curve fits were seen. The influence on the other photo- 
metric elements of adopting different limb-darkening coeffi- 
cients/laws are discussed below for the individual systems. 
In order to include uncertainties of the limb-darkening co- 
efficients in the uncertainties of the final photometric el- 
ements, random variations within +/- 0.05 have been in- 
troduced in the JKTEBOP Monte Carlo simulations men- 
tioned above. 
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In tables and text on photometric solutions we use the 
following symbols: i orbital inclination; r relative radius; 
k = r s /r p ; u linear limb darkening coefficient; y gravity 
darkening coefficient; J central surface brightness; L lumi- 
nosity; Q surface potential; T e g effective temperature. 

2.5. Effective temperatures and interstellar reddening 

Although significantly improved (semi)empirical tempera- 
ture scales and calibrations have become available during 
the last decade, the determination of effective tempera- 
tures of binary components remains a difficult task; see e.g. 
Popper PMTJ|) . Andersen (TTMTj) . and Clausen (|2"DU4) for 
general remarks. For comparison with theoretical models, 
luminosity calculations, and distance determinations, it is 
on the other hand important to diminish both systematic 
and random errors as much as possible. In the analyses of 
AD Boo, VZ Hya, and WZ Oph we have chosen to follow 
three independent approaches: 

First, as part of the light curve analyses we calculate the 
individual uvby indices for the components from the com- 
bined standard uvby indices (CVG08) and the luminosity 
ratios between the components in the four bands (based on 
mean geometric elements). Interstellar reddening is deter- 
mined from the intrinsic colour calibration by Olsen (1988), 
which has a precision of 0.009 mag, and the results are com- 
pared to those from other sources. Effective temperatures 
can then be obtained from a series of (semi) empirical cal- 
ibrations such as those by Alonso et al. (fT9~96h Ramirez 
& Melendez ([2005]) . and Holmberg et al. (|2"UU?j) . Also, as 
an important check on these calibrations, effective tempe- 
rature differences between the components can be derived 
from the empirical flux calibration by Popper (T98D|) and 
the y surface flux ratios. They are obtained to high preci- 
sion in the light curve analyses, since they depend mainly 
on the observed depths of primary and secondary eclipses. 
We have applied both the BC scale used by Popper (1980) 
and that by Flower (|1996p . which in the temperature range 
of AD Boo, VZHya, and WZOph, lead to nearly identi- 
cal temperature differences with errors of only about 25 K 
(excluding uncertainties in Popper's flux scale). 

Second, we attempt to obtain spectroscopic tempera- 
tures from the high-resolution and high-S/N spectra, either 
combined or disentangled, as part of the abundance analy- 
ses; see e.g. Santos et al. (|2004|) and Bruntt et al. (|2004|) 
for details. And thirdly, as mentioned in Sect. 12.21 tempera- 
tures are estimated by identifying the templates producing 
the best match to each component in the TODCOR analy- 
ses, and interpolating in log g or metallicity as necessary. As 
discussed below for the individual systems, the two latter 
methods however tend to be problematic. 

Finally, the photometric indices of the combined light 
of a system may provide a first estimate of the effective 
temperature of its "average" component, and if the com- 
ponents are nearly identical, as is the case for WZ Oph, a 
direct determination is possible. In Table[2]we present such 
temperatures derived from the combined uvby indices and 
2MASS0 JHK S observations. As seen, the temperatures 
from the different uvby and VJHK S calibrations agree well 
for WZ Oph. The spread for AD Boo and VZ Hya presu- 
mably reflects the fact that their combined flux distribu- 



tions are distorted compared to a single star, because the 
components differ by 400-600 K in effective temperature. 

2.6. Abundance analyses 

For meaningful comparisons between theoretical stellar 
models and empirical absolute dimensions — as determined 
from detached double-lined eclipsing binaries — the chemi- 
cal composition, and in particular the metal abundance, 
must be known for the binary. However, detailed abun- 
dance studies have only been done for a few systems. As 
part of the analyses of AD Boo, VZ Hya, and WZ Oph, we 
have used the new spectroscopic and photometric data to 
establish their metal abundance. Several methods were ap- 
plied/attempted: 

1. Spectroscopic analyses of observed or disentangled 
FEROS spectra 

2. Spectroscopic analyses of disentangled CfA spectra 

3. Template matching as part of the TODCOR analyses 

4. Photometric uvby-[Fe/H] calibrations 

2.6.1. FEROS spectra 

The FEROS spectra (Table [J) of WZ Oph were disentan- 
gled in order to extract the individual component spectra. 
We applied the disentangling method introduced by Simon 
& Sturm (1994J and a revised version of the corresponding 
original code developed by E. Sturm. It assumes a con- 
stant light level, so only spectra taken outside of eclipse 
can be included. We note that all three systems are in 
fact constant to within a few percent outside of eclipse. 
The orbital elements were fixed at the results based on 
the TODCOR analyses of the CfA spectra. For AD Boo, 
with only two FEROS spectra available, disentangling is 
not possible. Disentangling of the five VZ Hya spectra was 
attempted, but the quality of the resulting component spec- 
tra was not acceptable; more than five spectra and/or a 
better phase distribution seem to be needed. For both sy- 
stems the observed double- lined spectra were therefore used 
directly for the abundance analyses. 

The versatile semi-automatic IDL tool VW^Q develo- 
ped by Bruntt, was applied for the abundance analyses. 
We refer to Bruntt et al. (|2lK?2l I20W1 [211051) for further 
details. For our particular study, VWA has been extended 
to allow the analysis of double-lined spectra in addition to 
normal spectra from single stars — or disentangled binary 
component spectra. We have adopted the effective tempera- 
tures, surface gravities, and rotational velocities determined 
for the binary components in this study (Table l2"Tj) . but for 
comparison analyses for higher or lower temperatures were 
also done. Microturbulcncc velocities were derived from the 
calibration by Edvardsson et al. (|1993p , which has a scatter 
of about 0.3 kms -1 , but have been tuned, if needed, as part 
of the analyses. 

In each spectrum, VWA selects the least blended lines 
from atomic line lists extracted from the Vienna Atomic 
Line Database (VALD; Kupka et al. 11999j) . For each se- 
lected line, the synthetic spectrum is then calculated, and 
the input abundance is changed iteratively until the equiva- 
lent width of the observed and synthetic lines match. The 
synthetic spectra are generated with the SYNTH software 
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Table 2. Effective temperatures (K) of the 'average' components of AD Boo, VZHya, and WZOph. The Vb magnitudes 
and the (b—y)o and Co indices are based on the out-of-eclipse uvby standard indices from CVG08. The 2MASS observations 
(J, H, K a ) were obtained at phases 0.216 (AD Boo), 0.181 (VZ Hya), and 0.495 (WZ Oph), respectively. For WZ Oph, the 
V magnitude was calculated from the value outside eclipses and the y eclipse depth at phase 0.495. Ay is the adopted 
visual interstellar absorption. References are: A96, Alonso et al. (1996); H07, Holmberg et al. (|2007|) : RM05, Ramirez & 
Melendez ([2005]); M06, Masana et al. (|2006j) . The results from A96 are based on their uvby calibration, those from RM05 
on their uvby, (V — J), (V — H), and (V — K s ) calibrations (in that order); the calibration by M06 is for (V — K s ). 



System 


A v 


[Fc/H] 


(b - y) 


Co 


(V- J) 


(V-H) 


(V-K s ) 


A96 


H07 


RM05 


M06 


AD Boo 


0.107 


+0.10 


0.297 


0.426 


0.915 


1.141 


1.182 


6456 


6448 


6524/6370/6197/6296 


6333 


VZ Hya 


0.086 


-0.20 


0.280 


0.385 


0.817 


1.036 


1.096 


6511 


6517 


6528/6707/6463/6480 


6513 


WZ Oph 


0.141 


-0.27 


0.333 


0.369 


0.979 


1.176 


1.263 


6115 


6136 


6097/6190/6191/6176 


6210 



(Valenti & Piskunov ll996p . Atmosphere models were either 
interpolated from the grid of modified ATLAS9 models by 
Heiter et al. (|2002p or calculated from the MARCS code. 
We have applied an improved version of the original model 
by Gustafsson et al. (1975); see J0rgensen et al. (1992) for 
further details. 

For the binary components, we calculate all abundances 
relative to the Sun, line by line. We have used either an 
atlas spectrum of the solar flux (Wallace et al. 11998ft or 
a FEROS sky spectrum. In general, this approach gives 
better agreement between the abundances determined from 
the individual lines of a given element/ion. 

If we adopt the VALD line data and the solar photo- 
sphere abundances by Grevesse & Sauval (|1998p . we obtain 
for the Sun the Fe abundances listed in Table [3] for the so- 
lar atlas and FEROS spectra. Only lines with equivalent 
widths between 10 and 80 mA measured in both spectra 
were included. It is gratifying to see that the results from 
the FEROS spectrum agree very well with those from the 
atlas spectrum. For both spectra, Fel and Fell lines give 
slightly different [Fe/H] results. They agree a bit better 
for the MARCS models than for ATLAS9 models, and the 
MARCS models also reproduce the solar result slightly bet- 
ter. 

Although we will focus on Fe, we have also analysed 
other elements, and results from the solar atlas spectrum 
are given in Table [U For most of the ions with many mea- 
sured lines we obtain good agreement with Grevesse & 
Sauval (1998 ), but there are exceptions such as Si I, CrII, 
and Mnl. This may indicate errors in the VALD line data 
(a problem with Mn was found by Bruntt et al. l2004|) which 
should, however, be eliminated to a large extent in the bi- 
nary abundance analyses performed relative to the Sun. 
Comparison with the most recent solar photosphere abun- 
dances by Grevesse et al. (I2007p . which for many elements 
are derived using updated line data, 3D models and depar- 
tures from LTE, is less obvious, since we apply ID, LTE 
atmosphere codes. 



Table 3. Iron abundance ([Fe/H]) for the Sun deter- 
mined from Solar atlas spectrum (Wallance et al. [1998 ) 
and FEROS spectra, respectively. 131 Fe I and 20 Fe II lines, 
measured in both spectra, were used. Results from ATLAS9 
models and MARCS models are compared. 





Fel 


Fell 


ATLAS9: 






Solar atlas 


+0.11 ±0.10 


-0.02 ±0.07 


FEROS 


+0.11 ±0.10 


-0.04 ±0.07 


MARCS: 






Solar atlas 


+0.07 ±0.09 


±0.00 ±0.10 


FEROS 


±0.07 ±0.10 


-0.01 ±0.09 



luminosity ratios between the components, which cannot 
be determined from the disentangling procedure, were cal- 
culated by interpolation between the b and y luminosity 
ratios of the photometric solutions. The disentangled com- 
ponent spectra were then compared to the CfA library of 
synthetic spectra mentioned in Sect. 12.21 

In general we find that the short wavelength range cove- 
red (45 A) , and the correspondingly small number of lines 
available, limits the success of this approach. Most impor- 
tantly, it was difficult to break the correlation between 
metal abundance and effective temperature. Results for the 
individual systems will, however, be presented below. 

2.6.3. TODCOR template matching 

As mentioned in Sect. 12.21 rough estimates of the metal 
abundance for a system can be obtained directly from the 
observed composite spectra by running extensive grids of 
two-dimensional cross-correlations (TODCOR) using tem- 
plates of different metallicity. Results from this exercise will 
also be discussed below for the individual systems. 



2.6.2. Disentangled CfA spectra 

As seen in Tables [5j [TT] and [16l a large number of CfA 
spectra are available for all three systems. We have selected 
those collected outside eclipses with sufficient S/N ratios for 
disentangling. The actual numbers of spectra used are 113 
(AD Boo), 29 (VZHya), and 40 (WZOph). After careful 
normalization of the individual spectra, the disentangled 
component spectra were constructed by adopting the or- 
bital elements listed in the tables mentioned above. The 



2.6.4. uvby calibration 

Finally, we have derived metal abundances from the uvby 
indices for the individual binary components ( Table [2~lj) and 
the recent calibrations by Holmberg et al. (|20U7|) , which are 
based on a carefully selected sample of spectroscopic [Fc/H] 
results for several hundred stars. The calibration adopted 
for 0.30 < b — y < 0.46 has a dispersion of only 0.07 dex; 
that for bluer stars has a dispersion of 0.10 dex. Taking 
into account the uncertainties of the uvby indices and the 
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Table 4. Abundances for the Sun determined from Solar 
atlas spectrum (Wallace et al. I1998P and ATLAS9 models 
(Heiter et al. [2002 ) . N is the number of lines used; only ions 
with three or more measured lines are listed. GS1998 is the 
logarithmic clement abundances (H = 12.00) from Grevesse 
& Sauval (|T998)) . 
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Cr II 


0.16 


0.15 


22 


5.67 ±0.03 


Mn I 


0.19 


0.22 


39 


5.39 ±0.03 


Fe I 


0.11 


0.10 


238 


7.50 ±0.05 


Fe II 


-0.02 


0.08 


30 


7.50 ±0.05 


Co I 


0.03 


0.18 


49 


4.92 ± 0.04 


Ni I 


0.03 


0.10 


93 


6.25 ±0.04 


Cu I 


-0.17 


0.27 


4 


4.21 ±0.04 


Zn I 


0.09 


0.13 


4 


4.60 ±0.08 


Y II 


0.08 


0.17 


14 


2.24 ±0.03 


Zr I 


-0.23 


0.07 


3 


2.60 ±0.02 


Ba II 


0.23 


0.08 


3 


2.13 ±0.05 


Ce II 


0.11 


0.21 


8 


1.58 ±0.09 


Nd II 


0.13 


0.19 


5 


1.50 ±0.06 


Sm II 


0.05 


0.23 


3 


1.01 ±0.06 
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Fig. 3. Systematic errors in the raw TODCOR veloci- 
ties of AD Boo determined from simulations with synthetic 
binary spectra (filled circles: primary; open circles: sec- 
ondary). The differences are plotted both as function of or- 
bital phase (upper panel) and radial velocity (lower panel) , 
and have been applied to the measured velocities as correc- 
tions. Phase 0.0 corresponds to central primary eclipse. 



interstellar reddening, an overall precision of 0.10 to 0.15 
dex can be achieved per component. 

3. Spectroscopic and photometric analyses 

In this section we present in detail the spectroscopic and 
photometric analyses of AD Boo, VZ Hya, and WZ Oph, in- 
cluding abundances, and compare our results with those ob- 
tained in previous studies. Additional references to earlier 
work are given by CVG08. 

3.1. AD Boo 

AD Boo = BD±25°2800 = BV1 35 (P = 2 d 07) was disco- 
vered by Strohmcicr et al. (1956) to be an eclipsing binary, 
and Lacy (|1985j) found it to be double-lined. It has most 
recently been studied by Lacy ()1997a|) and Popper (1998a). 

3.1.1. Radial velocities and spectroscopic orbit for AD Boo 

Radial velocities for the components of AD Boo were de- 
termined by TODCOR analyses of 124 usable CfA spec- 
tra outside of eclipse, observed from January 1989 to June 
1997. Surface gravities (logc/) were fixed at 4.0 (primary) 



and 4.5 (secondary), and the heavy element abundance was 
assumed to be solar, the closest values in our grid of syn- 
thetic spectra to the actual values determined in Sect. 01 
We inferred temperatures of 6250±150 K for the primary of 
AD Boo (the more massive star) , and 5680 ± 300 K for the 
secondary, both significantly lower than the photometric 
values, because the templates have lower metallicity than 
measured for AD Boo; see Sects. 13.1.31 and |U The mea- 
sured usini values are 38 ± 2 kms -1 and 37 ± 5 kms -1 , 
respectively, with the secondary being considerably more 
uncertain due to its much weaker lines. The synchronous 
values are 39 kms -1 for the primary and 30 kms -1 for the 
secondary. Lacy (1997a) reported 38 ± 1 and 31 ± 1 kms -1 , 
whereas Popper (1998a) did not give any values. Templates 
from our library with parameters closest to the measured 
values were used to derive the radial velocities. 

In the case of AD Boo we found that the radial veloci- 
ty corrections for systematics described in Sect. 12. 2[ which 
are shown in Fig. [3j were quite significant, reaching up to 
10 kms -1 for the secondary star and up to 2 kms -1 for the 
primary. The effect on the absolute masses is large, ~6% for 
the primary and ~4% for the secondary. The final veloci- 
ties including corrections are listed in Table [A~T1 The light 
ratio derived from the spectra, including also corrections 
for systematics analogous to those derived for the veloci- 
ties, is £ s /ip — 0.37 ± 0.02 at the mean wavelength of our 
observations (5187 A). 
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Fig. 4. Spectroscopic orbital solution for AD Boo (solid 
line: primary; dashed line: secondary) and radial velocities 
(filled circles: primary; open circles: secondary). The dot- 
ted line (upper panel) represents the center-of-mass veloc- 
ity of the system. Phase 0.0 corresponds to central primary 
eclipse. 



The spectroscopic orbital solution for AD Boo is given 
in Table and the observations and computed orbit are 
shown graphically in Figure 2] along with the residuals for 
each star. For comparison, we have also disentangled 113 



Table 5. Spectroscopic orbital solution for AD Boo. T is 
the time of central primary eclipse. 



Parameter 



Value 



Adjusted quantities: 
K p (kms" 1 ) 
K s (kms" 1 ) 
7 (kms -1 ) 

Adopted quantities: 
P (days) 

T (HJD-2 400 000) 

e 

Derived quantities: 
M p sin 3 i (Mq) 
M a sin 3 i (Mq) 
a p sini (10 6 km) 
a s sini (10 km) 
asini (Rq) 

Other quantities pertaining to the fit: 

N ob3 

Time span (days) 
Op (kms -1 ) 
o s (kms -1 ) 



106.13 ±0.19 
124.11 ± 0.36 
-16.52 ± 0.13 

2.06880704 
49311.43169 
0.00 

1.4103 ±0.0089 
1.2061 ± 0.0056 
3.0193 ± 0.0053 
3.5307 ± 0.0103 
9.411 ±0.017 

124 
3074 
1.63 
3.15 



CfA spectra outside of eclipse and obtain velocity semi- 
amplitudes of K p = 105.95 and K s = 124.50, in good agree- 
ment with the results based on TODCOR velocities. 

Our velocity semiamplitudes agree well with those by 
Popper (11998afl t (K p K s ) = (105.85 ± 0.35, 124.4 ±0.4), 
whereas Lacy (|1997ap obtained (K P ,K S ) = (107.2 ± 0.5, 
124.6 ±0.6). 



3.1.2. Photometric elements for AD Boo 

The uvby light curves of AD Boo contain 652 observations 
in each band and were observed on 42 nights during five 
periods between March 1988 and March 1992 (CVG08). 
The accuracy per point is about 0.005 mag (vby) and 0.007 
mag (it), but comparison of the data from the five peri- 
ods reveals that at some phases, systematic differences in 
the light level of 0.01-0.03 mag exist, increasing from y to 
u. In particular, observations during the 1989 eclipses are 
systematically fainter. The effect is shown in Fig. [5] and in- 
dicates surface activity (spots) at a rather low level on one 
or both components, the cooler secondary being the most 
likely candidate. Zhai et al. (|19821) also noticed scatter in 
their B and V light curves, which were obtained from April 
to June 1981, but at a much higher level of about 0.075 mag 
compared to their observing error of about 0.02 mag. 

There is very little information on activity in the liter- 
ature, and AD Boo does not appear to be an X-ray source. 
The presence of chromospheric activity in terms of the 
Rossby numbeiQ was studied by Hall (|1994p . The secondary 
in AD Boo is perhaps cool enough to display significant ac- 
tivity (see above) despite having a mass as high as 1.21 
Mq, which would imply a relatively thin convective enve- 
lope. An estimate of the color index of B — V = 0.58 based 
on the temperature allows the convective turnover time to 
be inferred, and this together with the rotational period (as- 
suming synchronous rotation) leads to a Rossby number of 
approximately 0.20, which is in the range where other stars 
tend to show spot activity (see Hall 11994]) . Furthermore, 
our two high-resolution FEROS spectra taken in January 
1999 (see Table [lj show in fact a weak emission feature in 
the Ca II H and K lines at the position of the secondary 
component, supporting also some activity. Nothing is seen 
in the H Q line. 

EBOP and WD solutions, based on linear limb- 
darkening coefficients by van Hamme (1993) are presented 
in Tables [S] and [7J respectively, and O — C residuals of the 
v observations from the theoretical light curve computed 
for the photometric elements given in Table [5] are shown 
in Fig. O As seen, differences between the photometric el- 
ements obtained from the two models are small, and the 
mean values of r p + r s and k agree within their uncertain- 
ties. WINK solutions, which are not included, agree well 
with those from EBOP and WD. Slightly different orbital 
inclinations are obtained, however, for the four bands from 
all three codes. This could be due to either the systematic 
light level differences mentioned above (spots), which in- 
crease from y to u, to a small amount of (red) third light, 
or to a slightly incorrect theoretical wavelength dependence 
of the limb-darkening coefficients. 

Excluding the 1989 observations (see Fig [5]), we ob- 
tain the EBOP solutions listed in Table [8] The relative 



Defined as the ratio of the rotation period to the convective 
turnover time. 
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Table 6. Photometric solutions for AD Boo (all data) from 
the EBOP code adopting linear limb darkening coefficients 
by van Hamme (1993). The errors quoted for the free pa- 
rameters are the formal errors determined from the itera- 
tive least squares solution procedure. 





y 


b 


V 


it 


i (°) 


87.58 


87.69 


87.83 


87.89 


±12 


±10 


±10 


±15 


r p 


0.1711 


0.1715 


0.1709 


0.1713 




±5 


±4 


±4 


±7 


r s 


0.1291 


0.1292 


0.1289 


0.1296 


k 


0.754 


0.753 


0.754 


0.756 




±5 


±4 


±3 


±5 


r p + r s 


0.3002 


0.3007 


0.2998 


0.3009 


Up 


0.52 


0.61 


0.68 


0.65 


U s 


0.56 


0.65 


0.73 


0.75 


Vp 


0.34 


0.39 


0.45 


0.53 


Vs 


0.36 


0.42 


0.48 


0.57 


J s j Jp 


0.7314 


0.6888 


0.6440 


0.6775 


±33 


±30 


±31 


±48 


L>s/ Lp 


0.4079 


0.3828 


0.3566 


0.3693 


a (mag.) 


0.0064 


0.0063 


0.0067 


0.0102 



Table 7. Photometric solutions for AD Boo (all data) from 
the WD code adopting linear limb darkening coefficients by 
van Hamme (1993); see Table [51 Gravity darkening expo- 
nents of 0.33 and bolometric albedo coefficients of 0.5 were 
applied, as appropriate for convective envelopes. T c ff jP 
assumed to be 6575 K. 



was 





y 


b 


V 


it 


i n 


87.47 


87.78 


87.79 


87.82 


±5 


±3 


±3 


±6 




6.743 


6.711 


6.717 


6.702 




±11 


±8 


±9 


±14 


n s 


7.646 


7.696 


7.675 


7.636 




±16 


±10 


±9 


±18 


r v 


0.1705 


0.1713 


0.1711 


0.1716 


r s 


0.1303 


0.1293 


0.1297 


0.1305 


k 


0.764 


0.755 


0.758 


0.761 


r p + r s 


0.3008 


0.3006 


0.3008 


0.3021 


Teflf >a 


6078 


6052 


5999 


5946 


±4 


±4 


±4 


±4 




0.4205 


0.3868 


0.3615 


0.3761 


a (mag.) 


0.0064 


0.0062 


0.0066 


0.0101 



radii change by about 0.7% in opposite directions and their 
sum by about —0.6%, whereas the orbital inclination is 
unchanged in all four bands. Including 2% of third light 
changes i by about 0?3, which is the maximum difference 
between the four bands, but we have not found any direct 
support for even such a small amount of third light in any 
of the analyses. If we adopt linear limb-darkening coeffi- 
cients from Claret ([2000]) . which are 0.04-0.07 higher, the 
orbital inclination decreases by 0?16, and the relative ra- 
dius increases by about 0.0013 for the primary component 
but is unchanged for the secondary. If adjusted, the linear 
limb-darkening coefficients either become unrealistic and 
the orbital inclinations from the four bands disagree much 
more, or convergence fails. Finally, WD solutions including 
non- linear limb darkening (see Sect. 12 .4[) give photometric 
elements very close to those presented in Table [7J We are 



Table 8. Photometric solutions for AD Boo from the 
EBOP code, excluding observations from 1989. Linear 
limb darkening coefficients from van Hamme ( 1993) were 
adopted; see Table [5] 





V 




v 


u 


* ( ) 


87.59 


87.67 


87.83 


87.93 


±13 


±11 


±11 


±17 


r p 


0.1703 


0.1707 


0.1698 


0.1698 




±5 


±5 


±4 


±7 


r s 


0.1284 


0.1284 


0.1279 


0.1279 


k 


0.754 


0.752 


0.753 


0.753 




±6 


±4 


±3 


±6 


r p + r s 


0.2987 


0.2991 


0.2977 


0.2977 


J s j Jp 


0.7339 


0.6872 


0.6408 


0.6707 


±38 


±35 


±35 


±57 


L s 1 Lp 


0.4090 


0.3809 


0.3538 


0.3627 


a (mag.) 


0.0059 


0.0055 


0.0056 


0.0091 
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Fig. 5. (O—C) residuals of the AD Boo w-band observations 
from the theoretical light curve computed for the photomet- 
ric elements given in Table [6] 1989 observations are shown 
with square symbols. 



therefore unable to point at a clear explanation for the small 
i discrepancies, which are, however, of no consequence for 
the determination of the absolute properties. 

In order to assign realistic errors to the photometric e- 
lements, we have supplemented these analyses with 10 000 
JKTEBOP Monte Carlo simulations in each band; see 
Fig. [6] Furthermore, as shown in Fig. [7j we have com- 
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Table 9. Adopted photometric elements for AD Boo. The 
individual flux and luminosity ratios are based on the mean 
stellar and orbital parameters. 
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87? 73 ±0? 15 


r p 


0.1712 ±0.0015 


r„ 


0.1291 ± 0.0010 
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0.3003 ± 0.0020 
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0.754 ±0.006 
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Fig. 6. Best fitting parameter values for the 10 000 syn- 
thetic AD Boo z/-band light curves created for the Monte 
Carlo analysis. 
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Fig. 7. EBOP solutions for AD Boo for a range of fixed 
k values. Linear limb-darkening coefficients by van Hammc 
(1993) were adopted. The upper left panel shows normal- 
ized rms errors of the fit to the observations. Symbols are: 
cross y; diamond b; triangle v; square u. 



pared EBOP solutions for a range of fixed k = r s /r p va- 
lues near 0.75. The adopted photometric elements listed 
in Table [9] are the weighted mean values of the EBOP so- 
lutions adopting the linear limb-darkening coefficients by 
van Hamme. Errors are based on the Monte Carlo simula- 
tions and comparison between the uvby solutions, including 
those where the 1989 observations were excluded. We find 
that secondary eclipse is in fact total, whereas at phase 0.0 



J s I Jp 

L s J Lp 



V 

0.734 
±11 
0.410 

±6 



6 

0.690 
±11 
0.384 

±4 



0.642 
±10 
0.356 

±3 



0.678 
±11 
0.369 

±5 



about 61% of the y light from the primary component is 
eclipsed. The spectroscopic light ratio of £ s /£ p — 0.37±0.02 
at 5187A agrees quite well with that from the light curve 
solutions (0.40). 

For comparison, Lacy ( 1997aj), who applied the EBOP 
code to analyse the B and V normal points by Zhai et 
al. CEHE2), obtained i = 87?8 ± 0?2, r p = 0.1702 ± 0.0011, 
k = 0.750 ± 0.010, J s /J p = 0.719 ± 0.005 and L s /L p = 
0.403 ± 0.015 (both for the V band). 

3.1.3. Abundances for AD Boo 

Two high S/N FEROS spectra are available for abundance 
analyses of AD Boo; see Table [T] We have concentrated 
on the one with the highest S/N, which gives significantly 
more accurate abundances. Due to the relatively high ro- 
tational velocities and low luminosity ratio, the number of 
suitable unblended lines is far lower than for VZ Hya and 
WZ Oph, especially for the secondary component. ATLAS9 
models were used, and the effective temperatures, surface 
gravities and rotational velocities listed in Table [21] were 
adopted, together with microturbulcncc velocities of 1.5 
kms -1 for both components. The temperatures, which are 
derived from the uvby indices of the components, could not 
be constrained further from the spectroscopic analysis. 

As seen in Table [TU1 the [Fe/H] values from the Fel lines 
(both components) and the few Fell lines (primary) agree 
well, and identical results are obtained if MARCS models 
are used. Changing model temperatures by the uncertainty 
of the effective temperatures (±120 K) modifies [Fe/H] from 
the Fel lines by about ±0.10 dex, whereas the Fell result is 
changed by only ±0.04 dex. If 0.5 kms -1 higher microtur- 
bulence velocities are adopted, [Fe/H] decreases by about 
0.10 dex for both neutral and ionized lines. Taking these 
contributions into account, we obtain an average metallici- 
ty of [Fe/H] = ±0.07 ±0.14 for AD Boo. Abundances for Ca 
and Ni are based on very few lines; although higher than 
for Fe, they agree within the uncertainties. 

Metallicities have also been derived by comparing the 
disentangled CfA spectra to the extensive library of syn- 
thetic spectra especially prepared for the 45 A region cen- 
tered at 5187 A. Surface gravities were fixed as above, and 
synchronous rotation of the components was adopted. The 
best fits are obtained for [M/H] = ±0.13 and T off = 6450 K 
(primary), and [M/H] = ±0.30 and T eS = 6200 K (sec- 
ondary), but a strong correlation between metallicities and 
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Table 10. Abundances ([El./H]) for the primary and se- 
condary components of AD Boo. N is the number of lines 
used per ion. 





P 


rimary 




Secondary 




Ion 


[El./H] 


rms 


N 


[El./H] 


rms 


N 


Ca I 


+0.22 


0.06 


4 


+0.20 


0.23 


3 


Fe I 


+0.08 


0.15 


33 


+0.06 


0.17 


14 


Fe II 


+0.09 


0.21 


4 








Ni I 


+0.22 


0.11 


3 


+0.16 




1 



temperatures exists, and the determinations are therefore 
quite uncertain. We note that the temperature difference 
between the components becomes smaller than established 
from the light curve analyses. Forcing the metallicities of 
the components to be identical, we derive [M/H] = +0.25 
and temperatures of 6530 and 6160 K for the components, 
in good agreement with the photometric temperatures. 
Estimated uncertainties of metallicities and temperatures 
are (at least) 0.1 dex and 100 K. 

Grids of TODCOR analyses of the CfA spectra for a 
range of adopted template metallicities yield a metallicity 
estimate of [M/H] = +0.33, assuming temperatures of 6575 
and 6145 K, and forcing the results for the components to 
be identical. 

Finally, the calibrations by Holmberg et al. (|20U7|) and 
the de-reddened uvby indices (Table [H]) yield [Fe/H] = 
+0.28 ± 0.15 for the primary component ("blue" calibra- 
tion) and [Fe/H] = +0.10 ± 0.13 for the secondary. 

Within the fairly large errors, the results from the dif- 
ferent methods are in reasonable agreement. Giving higher 
weight to the analysis of the FEROS spectra, we adopt 
[Fe/H]= +0.10 ± 0.15 for AD Boo. This result should, 
however, be checked and preferably improved via analyses 
of disentangled component spectra based on a sufficiently 
large number of high-dispersion spectra. 

3.2. VZHya 

VZHya = HD72257 = HIP41834 (P = 2 d 90) was disco- 
vered by O'Connell ((T932) to be an eclipsing binary, and 
spectroscopic elements were obtained from 75 A/mm spec- 
tra by Struve (|1945p . The currently most reliable absolute 
dimensions of the components of VZ Hya are those given 
by Popper (|1980[) in his critical review on stellar masses. 

3.2.1. Radial velocities and spectroscopic orbit for VZ Hya 

Our TODCOR analysis of VZ Hya was based on 42 usable 
CfA spectra outside of eclipse, observed from January 1989 
to April 1991. Templates of solar metallicity, \ogg = 4.5 
(both components), and temperatures of 6750 K (primary) 
and 6250 K (secondary) were used to derive the radial velo- 
cities. As seen in Fig- [HI TODCOR corrections for systema- 
tics in the radial velocities are fairly small, typically under 1 
km s _1 for both components. They make very little differen- 
ce in the masses, changing them by only +0.4% (primary) 
and —0.1% (secondary), which is well within the errors. 
The scatter around the orbit is about the same with and 
without corrections applied. Interpolating between log g = 
4.0 and 4.5 grids to the exact gravities given in Table |2"T1 we 



6 2 



° -3 - 



> 

02 




-100 100 
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Fig. 8. Systematic errors in the raw TODCOR veloci- 
ties of VZ Hya, determined from simulations with synthetic 
binary spectra (filled circles: primary; open circles: sec- 
ondary). The differences are plotted both as function of or- 
bital phase (upper panel) and radial velocity (lower panel) , 
and have been applied to the measured velocities as correc- 
tions. Phase 0.0 corresponds to central primary eclipse. 



obtain temperatures of 6590 ± 100 K for the primary and 
6150 + 150 K for the secondary, both in excellent agreement 
with the photometric values; see Sect. HI The corresponding 
light ratio is t a Ji v = 0.53 ± 0.02 at the mean wavelength 
of our observations (5187 A). Finally, the grids of cross- 
correlations give best vsvcii values of 21 and 20 kms -1 , 
with errors of 2 kms -1 and 3 kms -1 , respectively; they are 
consistent with synchronous rotation. 

The spectroscopic orbital solution for VZHya is given 
in Table lll| and the observations and computed orbit are 
shown graphically in Figure [5] along with the residuals 
for each star. Our velocity semiamplitudes agree well with 
those by Popper (fT965|l . (K p ,K s ) = (94.3 + 2.3, 103.9 + 1.8) 
but are significantly more accurate. 

For comparison, we have also disentangled the 29 best 
CfA spectra outside of eclipse and obtain velocity semi- 
amplitudes of K p = 94.62 and K s = 105.87, in good agree- 
ment with the results based on TODCOR velocities. 



3.2.2. Photometric elements for VZ Hya 

The uvby light curves of VZHya contain 1180 observations 
in each band and were observed on 44 nights during four 
periods between February 1989 and April 1992 (CVG08). 
The average observational accuracy per point is about 4 
mmag (vby) and 6 mmag (u), and throughout all phases 
the points scatter at these levels, indicating that the com- 
ponents of VZ Hya are constant within the precision of our 
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Fig. 9. Spectroscopic orbital solution for VZ Hya (solid 
line: primary; dashed line: secondary) and radial velocities 
(filled circles: primary; open circles: secondary). The dot- 
ted line (upper panel) represents the center-of-mass veloc- 
ity of the system. Phase 0.0 corresponds to central primary 
eclipse. 



data. Furthermore, we see no signs of emission in CallH 
and K, or in Hq,, in the FEROS spectra. 

EBOP and WD solutions based on linear limb- 
darkening coefficients by van Hammc ( 1993 ) are shown in 



Table 11. Spectroscopic orbital solution for VZHya. T is 
the time of central primary eclipse. 



Parameter 



Value 



Adjusted quantities: 
K p (kms" 1 ) 
K s (kms" 1 ) 
7 (kms -1 ) 

Adopted quantities: 
P (days) 

T (HJD-2 400 000) 

e 

Derived quantities: 
M p sin 3 i (Mq) 
M a sin 3 i (Mq) 
a p sini (10 6 km) 
a s sini (10 km) 
asini (Rq) 

Other quantities pertaining to the fit: 

N ob3 

Time span (days) 
Op (kms -1 ) 
o s (kms -1 ) 



94.92 ± 0.19 
105.31 ±0.34 
-4.57 ± 0.14 

2.90430023 
48273.63450 
0.00 

1.2705 ± 0.0090 
1.1451 ±0.0061 
3.7907 ± 0.0079 
4.2057 ±0.0140 
11.489 ± 0.023 

42 
799 
1.04 
1.85 



Table 12. Photometric solutions for VZHya from the 
EBOP code adopting linear limb darkening coefficients by 
van Hamme (|1993p . The errors quoted for the free param- 
eters are the formal errors determined from the iterative 
least squares solution procedure. 





II 

y 


b 




u 


i (°) 


88.83 


88.91 


88.88 


88.89 


±6 


±5 


±5 


±8 


r p 


0.1142 


0.1144 


0.1144 


0.1140 




±4 


±3 


±3 


±5 


r s 


0.0969 


0.0967 


0.0969 


0.0970 


k 


0.849 


0.845 


0.847 


0.850 




±7 


±4 


±5 


±8 


r v ± r s 


0.2111 


0.2111 


0.2113 


0.2110 


Up 


0.52 


0.60 


0.67 


0.64 


U s 


0.54 


0.63 


0.71 


0.71 


Vp 


0.34 


0.39 


0.44 


0.52 


Vb 


0.35 


0.41 


0.47 


0.55 


J 3 1 Jp 


0.7983 


0.7665 


0.7306 


0.7649 


±20 


±16 


±18 


±29 


L s 1 Lp 


0.5701 


0.5402 


0.5143 


0.5364 


a (mag.) 


0.0046 


0.0040 


0.0044 


0.0068 



Table 13. Photometric solutions for VZHya from the 
WD code adopting linear limb darkening coefficients by 
van Hamme (| 1993[) ; see Table rT2] Gravity darkening expo- 
nents of 0.33 and bolometric albedo coefficients of 0.5 were 
adopted, as appropriate for convective envelopes. T ff jP was 
assumed to be 6650 K. 





V 


b 


V 


u 


i n 


88.90 


88.92 


88.95 


88.98 


±2 


±1 


±2 


±2 


flip 


9.649 


9.666 


9.639 


9.668 




±8 


±7 


±8 


±11 


n. 


10.349 


10.329 


10.355 


10.343 




±8 


±7 


±8 


±11 


r p 


0.1144 


0.1142 


0.1146 


0.1142 


r s 


0.0969 


0.0971 


0.0968 


0.0969 


k 


0.847 


0.850 


0.845 


0.849 


r v + r s 


0.2113 


0.2113 


0.2114 


0.2111 


T e ff, s 


6272 


6252 


6229 


6182 


±2 


±1 


±2 


±3 


Ls 1 Lp 


0.5686 


0.5464 


0.5125 


0.5347 


a (mag.) 


0.0046 


0.0039 


0.0043 


0.0068 



Tables fT2l and fl~3| respectively. As seen, the results from the 
different bands agree very well, and the two models lead to 
practically identical photometric elements. The same holds 
for WINK solutions not included here. If we adopt instead 
linear limb-darkening coefficients by Claret (2000), which 
are 0.05-0.08 higher, the orbital inclination decreases by 
0?2, and the relative radius is unchanged for the primary 
component but increases by 0.0009 for the secondary. If ad- 
justed, the linear limb-darkening coefficients by van Hamme 
are reproduced to within ±0.05 and the solutions are close 
to the results in Tables rT2] and [T3J WD solutions inclu- 
ding non- linear limb darkening (see Sect. 12 .4[) give photo- 
metric elements very close to those presented in Table rj31 
O — C residuals of the b observations from the theoretical 
light curve computed for the photometric elements given in 
Table [H] are shown in Fig. [TO] 
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Fig. 11. Best fitting parameter values for the 10 000 syn- 
thetic VZHya y light curves created for the Monte Carlo 
analysis. 



0.03 
0.02 



Fig. 10. (O — C) residuals of the VZHya 6-band obser- 
vations from the theoretical light curve computed for the 
photometric elements given in Table [T2"! 



As done for AD Boo, in order to assign realistic errors 
to the photometric elements, we have supplemented these 
analyses with 10 000 JKTEBOP Monte Carlo simulations in 
each band, and EBOP solutions for a range of fixed k values 
near 0.85; see Fig. [Til and Fig.[T2J The adopted photometric 
elements listed in Table [JJ] are the weighted mean values 
of the EBOP solutions adopting the linear limb-darkening 
coefficients by van Hamme. Errors are based on the Monte 
Carlo simulations and a comparison between the uvby solu- 
tions. Our results are significantly more accurate than ob- 
tained in previous studies of VZ Hya. At phase 0.0, about 
76% of the y light from the primary component is eclipsed, 
and at phase 0.5 about 98% of the y light of the secondary 
is blocked out, meaning that secondary eclipse is almost 
total. The spectroscopic light ratio of i s /i v — 0.53 ± 0.02 
at 5187A agrees quite well with that from the light curve 
solutions (0.56). 

Previous analyses based on the photoelectric UBV light 
curves by Walker (|1970p gave 3-4% larger relative radii for 
the primary component, whereas those for the secondary 
component agree well with our result (Walker [19701 Wood 
[T97T1 Cester et al. [T978| . Another set of photoelectric UBV 
light curves was published by Padalia & Srivastava (1975), 
who found primary eclipse to be an occultation (larger star 
in front) rather than a transit (smaller star in front). This 
picture of VZ Hya was soon shown by Popper (|1976p to be 
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Fig. 12. EBOP solutions for VZHya for a range of fixed 
k values. Linear limb-darkening coefficients by van Hamme 
(|1993p were adopted. The upper left panel shows normal- 
ized rms errors of the fit to the observations. Symbols are: 
cross y; diamond b; triangle v; square u. 



wrong but was nevertheless defended once more by Padalia 
(IT986I) . 



3.2.3. Abundances for VZ Hya 

Five FEROS spectra are available for the detailed abun- 
dance analysis of VZ Hya; see Table [1] As mentioned in 
Sect. I2.6[ disentangling of so few spectra was unsuccess- 
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Table 14. Adopted photometric elements for VZ Hya. The 
individual flux and luminosity ratios are based on the mean 
stellar and orbital parameters. 



i 

r p 



r v ± r s 
k 



88? 88 ±0? 09 
0.1143 ±0.0004 
0.0968 ± 0.0006 
0.2111 ±0.0004 
0.847 ± 0.007 



Js I ' Jjj 

L s /Lp 



y 


b 


V 


u 


0.799 


0.766 


0.730 


0.764 


±12 


±12 


±11 


±12 


0.569 


0.543 


0.515 


0.532 


±9 


±6 


±6 


±9 



Table 15. Abundances ([El./H]) for the primary and se- 
condary components of VZ Hya. N is the number of lines 
used per ion. 





P 


rimary 




Secondary 




Ion 


[El./H] 


rms 


N 


[El./H] 


rms 


N 


Si I 


-0.14 


0.17 


7 


-0.18 


0.16 


12 


Ca I 


-0.13 


0.15 


8 


-0.17 


0.16 


7 


Sc II 


-0.36 


0.14 


2 


-0.33 


0.07 


3 


Ti I 


-0.25 




1 


-0.21 


0.17 


4 


Ti II 


-0.48 


0.43 


3 


-0.30 


0.09 


2 


Cr I 


-0.23 


0.16 


6 


-0.24 


0.11 


4 


Cr II 


-0.29 


0.03 


4 


-0.43 


0.06 


4 


Fe I 


-0.21 


0.16 


91 


-0.19 


0.18 


69 


Fe II 


-0.32 


0.18 


12 


-0.14 


0.26 


11 


Ni I 


-0.20 


0.23 


9 


-0.07 


0.27 


16 



ful, and we have instead used the three most suitable ob- 
served spectra and analysed them as double-lined. The re- 
sults from the VWA analysis are presented in Table [TS] 
ATLAS9 models were used, and the final effective tempe- 
ratures, surface gravities and rotational velocities listed in 
Table |2"T1 were adopted, together with microturbulence ve- 
locities of 2.0 and 1.6 kms -1 , respectively, calculated from 
the calibration by Edvardsson et al. (|1993|) . For this choice 
of parameters there is no dependence on equivalent width 
and excitation potential for abundances determined from 
individual lines, which on the other hand appears if tem- 
peratures are changed by about 150 K. 

As seen, the [Fe/H] values from the many Fel lines 
agree well, whereas those from the much fewer Fell lines 
are lower (primary) and higher (secondary). Similar re- 
sults are obtained if MARCS models are used, and we 
have not been able to identify the cause of the discrepan- 
cies. Changing model temperatures by the uncertainty of 
the effective temperatures (±150 K) modifies [Fe/H] from 
the Fel lines by ±0.10 dex, whereas the Fell results are 
changed by only ±0.04 dex. If 0.5 kms -1 higher microtur- 
bulence velocities are adopted, [Fe/H] decreases by about 
0.06 dex for both neutral and ionized lines. Taking these 
contributions into account, we obtain a weighted average 
of [Fe/H] = -0.20 ±0.12 for VZHya. 

Excluding ions with few measurements, we find that 
within the uncertainties, identical abundances are derived 



from CrI and Nil lines. The a-elements Si and Ca seem 
slightly enhanced, but the difference is barely significant. 

For the metallicities derived by comparing the disen- 
tangled CfA spectra with the library of synthetic spectra, 
surface gravities were fixed as above, and synchronous ro- 
tation of the components was adopted. The best fits are 
obtained for [M/H] = ±0.09 and T cB = 6630 K (primary), 
and [M/H] = -0.23 and T off = 6130 K (secondary), but as 
for AD Boo, a strong correlation between metallicities and 
temperatures exist, and the determinations are therefore 
quite uncertain. Forcing the metallicities of the components 
to be identical, we derive [M/H] = —0.15, and temperatures 
of 6595 and 6185 K for the components, in good agreement 
with the FEROS result and the photometric temperatures. 
Estimated uncertainties of these CfA metallicities and tem- 
peratures are (at least) 0.1 dex and 100 K, respectively. 

Comparing TODCOR analyses of the observed CfA 
spectra for a range of adopted template metallicities yields 
a metallicity estimate of [M/H] = ±0.02, assuming temper- 
atures of 6650 and 6300 K and forcing the metal abundance 
of the components to be identical. 

Finally, the calibrations by Holmberg et al. (|2007p and 
the de-reddened uvby indices in Table [5T] yield [Fe/H] = 
±0.05±0. 15 for the primary component ("blue" calibration) 
and [Fe/H] = —0.08 ± 0.12 for the secondary, marginally 
in agreement with the FEROS result, which we adopt for 
VZ Hya. 

3.3. WZOph 

WZOph = HD154676 = HIP83719 (P = 4 d 18) was discov- 
ered by J. H. Metcalf to be an eclipsing binary (Pickering 
I1917p . and Sanford (|1937|) obtained the first radial ve- 
locities. The currently best modern analysis of WZ Oph, 
leading to masses and radii accurate to 3-4%, was done 
by Popper <| 1965|) from photoelectric light curves and 20 
A/mm spectra. 

For the spectroscopic and photometric analyses of 
WZ Oph we have adopted the well established ephemeris 
from CVG08, where the zero point has been shifted by half 
a period compared to that by Popper (|1965p . We note that 
the primary component eclipsed at phase 0.0 according to 
our adopted ephemeris, and corresponding to the slightly 
deeper minimum for the uvby light curves, is given by 
Popper (1965) as component 2. Popper found that his com- 
ponent 2, i.e. our primary component, has slightly stronger 
spectral lines. This is consistent with our measurements of 
equivalent widths in the disentangled FEROS spectra for 
lines stronger than about 100 mA. 

3.3.1. Radial velocities and spectroscopic orbit for WZOph 

Radial velocities for the components of WZ Oph were de- 
termined with TODCOR on the basis of 40 usable CfA 
spectra outside of eclipse, observed from February 1989 to 
September 1991. Identical templates of solar composition 
with log<? = 4.0 and a temperature of 6250 K were adopted 
for the two stars. The TODCOR corrections for systematics 
shown in Fig. [13] are less than 0.5 kms -1 at all phases for 
both components. This makes a difference of only ±0.3% 
in the minimum masses, about half their uncertainties, and 
the scatter of the orbit does not change much with correc- 
tions, wsini values of 18 and 17 kms -1 are obtained for the 
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Fig. 13. Systematic errors in the raw TODCOR veloc- 
ities of WZ Oph, determined from simulations with syn- 
thetic binary spectra (filled circles: primary; open circles: 
secondary). The differences are plotted both as function 
of orbital phase (upper panel) and radial velocity (lower 
panel), and have been applied to the measured velocities 
as corrections. Phase 0.0 corresponds to central primary 
eclipse. 
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Fig. 14. Spectroscopic orbital solution for WZ Oph (solid 
line: primary; dashed line: secondary) and radial velocities 
(filled circles: primary; open circles: secondary). The dot- 
ted line (upper panel) represents the center-of-mass veloc- 
ity of the system. Phase 0.0 corresponds to central primary 
eclipse. 



3.3.2. Photometric elements for WZOph 

The uvby light curves of WZ Oph contain 697 observations 
in each band and were observed on 41 nights during six 



primary and secondary components, respectively, with esti- 
mated uncertainties of 1 kms -1 . They are consistent with 
synchronous rotation. 

Our best spectroscopic estimates of the effective temper- 
atures of the components by interpolating to a log g value 
close to that listed in Table |2"T1 are 6260 K for the primary 
and 6340 K for the secondary, assuming solar metallicity. 
Interpolating to the metallicity of [M/H] = —0.27 (see be- 
low) gives 6050 K and 6120 K. The lower values for the pri- 
mary compared to the secondary disagree formally with the 
trend from photometry, which gives a slightly hotter pri- 
mary based on the slightly deeper primary eclipse. However, 
the difference is well within the spectroscopic error estimate 
of 150 K. The light ratio is £ s /i p = 0.98 ± 0.02 at the mean 
wavelength of our observations (5187 A). 

The spectroscopic orbital solution for WZ Oph is given 
in Table [16] and the observations and computed orbit are 
shown graphically in Figure [H] along with the residuals for 
each star. For comparison, we have also disentangled the 
40 CfA spectra outside of eclipse and obtain velocity semi- 
amplitudes of K p — 88.62 and K s = 88.92, in good agree- 
ment with the results based on the TODCOR velocities. 

Our velocity semiamplitudes are larger than those by 
Popper (T19651) . {K p , K s ) = (85.9 ± 1.8, 87.0 ± 2.0), and as 
in the case of VZ Hya they are of much higher precision. 



Table 16. Spectroscopic orbital solution for WZ Oph. T 
is the time of central primary eclipse. 



Parameter 


Value 


Adjusted quantities: 




K p (kms" 1 ) 


88.77 ±0.19 


K s (kms -1 ) 


89.26 ± 0.24 


7 (kms -1 ) 


-28.36 ±0.13 


Adopted quantities: 




P (days) 


4.18350681 


T (HJD-2 400 000) 


50353.78331 


e 


0.00 


Derived quantities: 




M p sin 3 i (Mq) 


1.2262 ± 0.0073 


M s sin 3 i (Mq) 


1.2194 ± 0.0063 


a p sini (10 6 km) 


5.106 ±0.011 


a s sin i (10 6 km) 


5.135 ±0.014 


asini (Rq) 


14.714 ±0.026 


Other quantities pertaining to the fit: 




N obs 


40 


Time span (days) 


918 


(jp (kms -1 ) 


1.04 


a s (kms - ) 


1.33 
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Table 17. Photometric solutions for WZOph from the 
EBOP code adopting linear limb darkening coefficients by 
van Hammc (1993). The errors quoted for the free param- 
eters are the formal errors determined from the iterative 
least squares solution procedure. 





y 


b 


V 


u 


i n 


89.08 


89.09 


89.06 


89.11 


1 O 


1 O 

ztz 


in 

±Z 


_|_ C 

±0 


r p 


0.0956 


0.0950 


0.0951 


0.0949 




±8 


±7 


±7 


±11 


r s 


0.0963 


0.0963 


0.0960 


0.0975 


k 


1.008 


1.014 


1.009 


1.027 




±16 


±14 


±14 


±22 


r p + r s 


0.1919 


0.1913 


0.1911 


0.1924 


Up — Us 


0.56 


0.65 


0.73 


0.74 


Vp = Vs 


0.36 


0.42 


0.48 


0.57 


J s 1 Jp 


0.9878 


0.9815 


0.9831 


0.9668 


±38 


±41 


±46 


±73 


L s f L p 


1.0035 


1.0102 


1.0001 


1.0209 


a (mag.) 


0.0053 


0.0057 


0.0065 


0.0104 



Table 18. Photometric solutions for WZ Oph from the 
WD code adopting linear limb darkening coefficients by 
van Hamme (|1993p ; see Tabic |T7J Gravity darkening expo- 
nents of 0.33 and bolometric albedo coefficients of 0.5 were 
adopted, as appropriate for convective envelopes. T g^ p was 
assumed to be 6165 K. 





y 


b 


V 


it 


i n 


89.11 


89.11 


89.08 


89.14 


±1 


±1 


±2 


±3 


alp 


11.454 


11.510 


11.554 


11.444 




±11 


±15 


±14 


±32 


Q, 


11.334 


11.299 


11.310 


11.276 




±12 


±14 


±17 


±30 


r p 


0.0957 


0.0952 


0.0948 


0.0958 


r s 


0.0963 


0.0967 


0.0966 


0.0969 


k 


1.007 


1.016 


1.019 


1.012 


r p + r 3 


0.1920 


0.1919 


0.1914 


0.1927 


T e ft,s 


6149 


6141 


6141 


6121 


±2 


±2 


±3 


±5 




1.0023 


1.0126 


1.0180 


0.9936 


a (mag.) 


0.0053 


0.0057 


0.0065 


0.0104 



periods between March 1991 and March 1997 (CVG08). 
The average observational accuracy per point is about 4 
mmag (vby) and 6 mmag (u), but throughout all phases 
the points scatter by 1-4 mmag more than this, highest in 
u. At a given phase, the observations from different seasons 
do not differ systematically, and it is thus unclear if the ex- 
tra scatter in the light curves is due to slight variability of 
one or both stars (activity), or perhaps to an unknown in- 
strumental/observational effect not seen in other systems. 
The Rossby number inferred for the mean component of 
WZ Oph does not clearly place it in the group of stars that 
show variability due to spots (Hall 1994) , and is thus incon- 
clusive. On the other hand, the high S/N FEROS spectrum 
taken at phase 0.33 (see Tabled]) indicates very weak emis- 
sion from both components in the Ca 11 H and K lines, but 
not in H a . This suggests WZ Oph is in fact mildly active. 

EBOP and WD solutions, based on linear limb- 
darkening coefficients by van Hammc (1993) are presented 



WZ Oph u 



r 


+ y~ ,4. \- - 


* 


- 










-0.03 -0.02 


-0.01 0.00 0.01 
Phase 


0.02 0.03 




++ 


i- ++ V + + + ++ 




+ + ^++ + + 1 






+ +++ _I 


0.47 0.48 


0.49 0.50 0.51 
Phase 


0.52 0.53 


- 

; :A f| tij 


% ^+ * * 




+ 






"1 * + 


+ 


0.1 


0.2 0.3 
Phase 


0.4 






t * t* 


P^h^i 




++ ++ + 




* + _ 


0.6 


0.7 


0.8 


0.9 



Fig. 15. (O — C) residuals of the WZOph y-band obser- 
vations from the theoretical light curve computed for the 
photometric elements given in Table [T71 

in Tables [17] and [TH] respectively. As seen, the results from 
the different bands agree well, and the two models lead to 
nearly identical photometric elements. WINK solutions not 
included here lead to changes in the relative radii of —0.0008 
(primary) and ±0.0008 (secondary). Formally, the less mas- 
sive secondary component is slightly larger and cooler, but 
more luminous, than the primary for all three models. In re- 
ality the two stars are, however, very nearly identical. If we 
adopt linear limb-darkening coefficients by Claret (2000), 
which are 0.10 higher than those by van Hamme, the or- 
bital inclination decreases by 0?1, and the relative radius 
increases by about 0.0004 for the primary component but 
is unchanged for the secondary. If adjusted, the vby linear 
limb-darkening coefficients by van Hamme are reproduced 
to within ±0.05 and the solutions are close to the results 
in Table 1171 The u coefficients, however, become unrealis- 
tically low (0.60), and i changes by 0?2. WD solutions in- 
cluding non- linear limb darkening (see Sect. 12 .4[) give pho- 
tometric elements very close to those presented in Table [TH] 
for the uvb bands, but the u solution now disagrees in hav- 
ing k w 1.05. O — C residuals of the y observations from 
the theoretical light curve computed for the photometric 
elements given in Table [171 are shown in Fig. [T5l 

Realistic errors have been assigned to the photometric 
elements by performing 10 000 JKTEBOP Monte Carlo 
simulations in each band, and EBOP solutions for a range of 
fixed k values near 1.0; see Fig. [TBI and Fig.QTJ The adopted 
photometric elements listed in Table [T5] are the weighted 
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| 750 

% 700 
1 650 
600 
550 




0.96 0.98 1.00 1.02 1.04 1.06 




1092 0.094 0.096 0.C 



Table 19. Adopted photometric elements for WZOph. 
The individual flux and luminosity ratios are based on the 
mean stellar and orbital parameters. 



i 


89? 08 ±0? 10 


r p 


0.0952 ± 0.0008 


r„ 


0.0964 ± 0.0008 


r p ± r s 


0.1916 ±0.0004 


k 


1.013 ±0.015 



89.20 
89.15 
89.10 
89.05 
89.00 




1.10 
1.05 



/ 



0.96 0.98 1.00 1.02 1.04 1.06 



0.96 0.98 1.00 1.02 1.04 1.06 



Fig. 16. Best fitting parameter values for the 10000 syn- 
thetic WZ Oph y light curves created for the Monte Carlo 
analysis. 



1.000 
0.995 
0.990 
0.985 
0.980 
0.975 
0.970 





0.96 0.98 1.00 1.02 1.04 1.06 
k = r_s/r_p 



0.96 0.98 1 .00 1 .02 1 .04 1 .06 
k = r_s/r_p 



0.098 
0.097 
0.096 
0.095 
0.094 
0.093 




0.96 0.98 1.00 1.02 1.04 1.06 
k = r_s/r_p 




).96 0.98 1.00 1.02 1.04 1.06 
k = r_s/r_p 



o 1.000 
» 0.990 



0.970 
0.960 





0.96 0.98 1.00 1.02 1.04 1.06 
k e r_s/r_p 



0.96 0.98 1 .00 1 .02 1 .04 1 .06 
k = r_s/r_p 



Fig. 17. EBOP solutions for WZ Oph for a range of fixed 
k values. Linear limb-darkening coefficients by van Hammc 
( 1993]) were adopted. The upper left panel shows normal- 
ized rms errors of the fit to the observations. Symbols are: 
cross y; diamond 6; triangle v; square u. 



mean values of the EBOP solutions adopting the linear 
limb-darkening coefficients by van Hamme. Errors are based 
on the Monte Carlo simulations and comparison between 
the uvby solutions. Adopting the spectroscopic light ratio 
of £ s /l p = 0.98 ± 0.02 at 5187 A would lead to k = 0.996 
(formal error ±0.011), r p = 0.0960, and r s = 0.0956. We 



y 

Js/Jp 0.987 
±15 

L a /L„ 1.013 
±30 



6 

0.981 
±16 
1.007 

±27 



0.982 0.972 

±16 ±17 

1.008 0.997 

±26 ±37 



find that about 92% of the light of the components is shad- 
owed out at both central eclipses. 

For comparison, Popper (|1965|) obtained slightly smaller 
relative radii (r p = 0.0945 ± 0.002, r s = 0.092 ± 0.001), and 
an orbital inclination of i = 89?4 ± 0?2 (all uncertainties 
given by Popper are estimated). From a WINK re-analysis 
of Popper's photoelectric light curve, Cester et al. {1978) 
derived r p = 0.096 ± 0.001, k = 0.971 ± 0.017, and i = 
89?0±0?1. 



3.3.3. Abundances for WZOph 

14 FEROS spectra are available of WZOph (see Table [J), 
which allowed us to obtain disentangled spectra of high 
quality for the abundance studies. Results from the VWA 
analysis are presented in Table [20] ATLAS9 models were 
adopted, and the effective temperatures, surface gravities 
and rotational velocities listed in Table [21] were used, 
together with microturbulence velocities of 1.7 kms -1 
for both components, calculated from the calibration by 
Edvardsson et al. (1993]). This choice of parameters is sup- 
ported by the spectroscopic analyses, since we see no de- 
pendence on equivalent width and excitation potential for 
abundances determined from individual lines. The precision 
of the temperatures (±100 K) could not be improved from 
the abundance analysis. 

As seen, a robust [Fe/H] is obtained with identical re- 
sults for Fel and Fell lines from both components, also if 
MARCS models are used. Changing model temperatures 
by ±100 K modifies [Fe/H] from the Fel lines by ±0.06 
dex, whereas the Fell results are nearly unchanged. If 
0.5 kms -1 higher microturbulence velocities are adopted, 
[Fe/H] decreases by about 0.04 dex for both neutral and 
ionized lines. Taking these contributions into account, we 
find [Fc/H] = -0.27 ± 0.07 for WZ Oph. 

Excluding ions with only a few lines measured, we find 
that within the uncertainties identical abundances are de- 
rived for Cr and Ni, noting that CrII deviates, as seen 
for the Sun (Table [4]), although our abundances are rela- 
tive so that the effect should cancel out. The a-elemcnts 
Si and Ca seem to be enhanced by about 0.1 dex, whereas 
Na and Ti (with only few lines measured) appear normal. 
The Si result may perhaps be questioned, as in the case 
of CrII above, since Si I gave a low abundance for the Sun. 
For Mn, which together with Al is expected to anticorrelate 
with a-elements, we actually derive a rather low abundance 
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Table 20. Abundances ([El./H]) for the primary and se- 
condary components of WZ Oph. N is the number of lines 
used per ion. 





P 


rimary 




Secondary 




Ion 


[El./H] 


rms 


N 


[El. /HI 

V I J 


rms 


N 


Nal 


-0.32 


0.05 


4 


—0.35 


0.05 


3 


Si I 


-0.19 


0.12 


20 


—0.18 


0.09 


16 


Ca I 


-0.21 


0.13 


11 


-0.36 


0.18 


9 


Sc II 


-0.28 


0.14 


5 


-0.30 


0.07 


5 


Ti I 


—0.32 


0.23 


3 


-0.26 


0.20 


3 


Ti II 


-0.12 


0.03 


4 


-0.17 


0.10 


4 


V I 


+0.03 


0.25 


10 


+0.06 


0.16 


3 


Cr I 


-0.36 


0.25 


11 


-0.39 


0.25 


10 


Cr II 


-0.17 


0.17 


8 


-0.14 


0.13 


6 


Mn I 


-0.54 


0.27 


9 


-0.64 


0.22 


5 


Fe I 


-0.27 


0.12 


60 


-0.28 


0.16 


74 


Fe II 


-0.28 


0.04 


9 


-0.27 


0.08 


9 


Co I 


+0.08 


0.13 


4 


-0.18 


0.31 


2 


Ni I 


-0.25 


0.23 


13 


-0.25 


0.16 


15 



of —0.54+0.13, taking into account temperature and micro- 
turbulence uncertainties as above. Unfortunately, problems 
with the Mnl line were also seen for the Sun. We tend, how- 
ever, to conclude, that a slight a-enhancement of about 0.1 
dex may be present. Finally, we derive an abundance for V 
close to solar. 

As for AD Boo and VZ Hya, metallicities have also been 
derived from the disentangled CfA spectra. Surface grav- 
ities were fixed as above, and synchronous rotation of 
the components was adopted. The best fits are obtained 
for [M/H]= -0.14 and T cff = 6280 K (primary), and 
[M/H] = -0.29 and T off = 6250 K (secondary), but as in 
the other two binaries there is a strong correlation between 
metallicity and temperature, and the determinations are 
therefore quite uncertain. Forcing the metallicities of the 
components to be identical, we derive [M/H] = —0.23 and 
temperatures of 6230 and 6310 K for the components, in 
good agreement with the FEROS abundance result and in 
reasonable agreement with the photometric temperatures. 
Estimated uncertainties of metallicities and temperatures 
from this method are (at least) 0.1 dex and 100 K, respec- 
tively. 

TODCOR analyses of the observed (composite) CfA 
spectra for a range of fixed template metallicities yield an 
estimate of [M/H] = —0.27, assuming temperatures of 6165 
and 6115 K and constraining the results for the components 
to be identical. 

Finally, the calibration by Holmberg et al. (|2007D and 
the de-reddened uvby indices yield [Fe/H] = —0.17 ± 0.11 
and —0.24 ± 0.12 for the primary and secondary compo- 
nents, respectively, also in good agreement with the FEROS 
result, which we adopt for WZ Oph. 

4. Absolute dimensions 

Absolute dimensions for AD Boo, VZ Hya, and WZ Oph are 
presented in Table [HI as calculated from the new spectro- 
scopic and photometric elements given in Tables [5j O HU 
[HQ1 and[H] 

As seen, masses and radii precise to 0.5-0.7% and 0.4- 
0.9%, respectively, have been established for the binary 



components. For AD Boo, the masses agree well with those 
determined by Popper (|1998ap but are significantly more 
accurate, whereas those by Lacy (1997a) are about 2% 
higher. Both studies report radii in agreement with ours. 
Popper (1965] I1980P reached masses and radii accurate to 
about 3% for VZHya, and within this uncertainty his re- 
sults agree with ours. The less massive secondary compo- 
nent of WZ Oph appears slightly larger and more luminous 
but cooler than the primary. However, within the uncer- 
tainties the two components can be considered identical. 
On average, we obtain masses and radii for WZ Oph which 
are about 8% and 5% larger, respectively, than those given 
by Popper (|1965j) . and they are far more accurate. 

The uvby indices for the components of AD Boo, 
VZ Hya, and WZ Oph, included in Tabled]] were calculated 
from the the combined indices of the systems at phase 0.25 
(CVG08) and the luminosity ratios between their compo- 
nents (Tables El [H HU). The E(b - y) interstellar redden- 
ings, also given in Table [2T1 were determined from uvby (3 
standard photometry and other sources, as described in 
more detail in Sect. 14.11 

The adopted effective temperatures listed in Table [2T] 
were calculated from the calibration by Holmberg et al. 
(2007) assuming the final [Fe/H] abundances. The uncer- 
tainties (100-150 K) include those of the uvby indices, 
E(b — y), [Fe/H], and the calibration itself. Within errors, 
temperatures from the (b — y), c\ calibration by Alonso et 
al. (|1996p and the b — y calibration by Ramirez & Melendez 
(2005J agree, but for AD Boo we note that both calibra- 
tions lead to a temperature difference between the com- 
ponents that is about 120 K higher. A similar trend was 
noticed by Holmberg et al. (2007). The empirical flux scale 
by Popper (|1980p and the y flux ratios between the com- 
ponents (Tables [H [T31 yield well established tempera- 
ture differences of 450 K (AD Boo), 335 K (VZHya), and 
20 K (WZ Oph), supporting the temperature calibration we 
have chosen. Details on temperatures derived as part of the 
TODCOR analyses of the observed CfA spectra (assuming 
the final metallicities) and the abundance analyses of the 
disentangled CfA spectra are given in Sect. [31 In general, 
the results agree well with those included in Table [HI ex- 
cept for WZ Oph where slightly higher temperatures were 
obtained from the analysis of the disentangled spectra. 

The distances to the binaries were calculated from the 
"classical" relation, 

V - My - (V - Ay) - M bolQ + 5 log(i?/i? Q ) (1) 
+ 101og(T cff /T cffQ ) + BC7 

adopting the solar values and bolometric corrections given 
in Table [U and A v /E(b - y) = 4.27. As seen, the dis- 
tances have been established to 5-6 %, accounting for all 
error sources and including the use of other BC scales 
(e g Code et al. IT9761 Bessell et al. IT9981 Girardi et al. 
2002). For AD Boo and VZHya, the empirical K surface 
brightness - T c g relation by Kervella et al. (2004) leads to 
nearly identical and perhaps even more precise distances; 
see Southworth et al. (2005) for details. K photometry out- 
side eclipses is unfortunately not available for WZ Oph. 
VZ Hya and WZ Oph are included in the Hipparcos ca- 
talog (ESA [TMTj) . which gives parallaxes of 5.03 ± 1.34 
mas and 7.99 ± 1.37 mas, respectively, corresponding to 
distances of 199 + 53 and 125 + 21 pc. Matching the shorter 
Hipparcos distance for WZ Oph would require lower tern- 
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Table 21. Astrophysical data for AD Boo, VZHya, and WZOph. T e//Q = 5780 K, M bolQ = 4.74, and bolometric 
corrections (BC) by Flower (|1996[) have been applied. v sync is the equatorial velocity for synchronous rotation. 





AD Boo 
Primary Secondary 


VZ Hya 
Primary Secondary 


WZ Oph 
Primary Secondary 


Ar\G*~»lii^~£a m on ci nn u ■ 
I\. USUI U. Lc LLllllfc/llolUllo . 














1V1 1 1V1 (.) 


1.414 ±0.009 


1.209 ± 0.006 


1.271 ±0.009 


1.146 ±0.006 


1.227 ±0.007 


1.220 ±0.006 


R/R,~, 


1.612 ± 0.014 


1.216 ± 0.010 


1.314 ±0.005 


1.112 ±0.007 


1.401 ± 0.012 


1.419 ±0.012 




4.173 ±0.008 


4.351 ± 0.007 


4.305 ± 0.005 


4.405 ± 0.006 


4.234 ± 0.008 


4.221 ±0.008 


v sin i (km s ) 


38 ±2 


37 ±5 


21 ±2 


20 ±3 


18 ± 1 


17 ±1 


Vsync (kHIS ) 


39.4 ±0.4 


29.7 ±0.2 


22.9 ±0.1 


19.4 ±0.1 


16.9 ±0.1 


17.2 ±0.1 


Photometric data: 
















9.752 ± 0.008 


10.721 ± 0.013 


9.442 ± 0.008 


10.054 ±0.012 


9.856 ± 0.017 


9.841 ±0.017 


v> - v)\ 


0.302 ± 0.006 


0.371 ± 0.007 


0.281 ± 0.005 


0.333 ± 0.005 


0.362 ± 0.005 


0.370 ± 0.005 


mit 


0.183 ±0.011 


0.198 ±0.011 


0.165 ± 0.008 


0.169 ±0.008 


0.149 ±0.008 


0.141 ±0.008 




0.464 ± 0.015 


0.340 ± 0.016 


0.420 ± 0.009 


0.329 ±0.010 


0.369 ±0.011 


0.382 ±0.012 


771/1, „ \ 

£,(6 - y) 


0.025 : 


t 0.015 


0.020 ± 0.020 


0.033 ±0.012 




6575 ± 120 


6145 ± 120 


6645 ± 150 


6290 ± 150 


6165 ± 100 


6115 ±100 


M bo j 


3.14 ±0.08 


4.05 ± 0.09 


6645 ± 150 


6290 ± 150 


3.73 ± 0.07 


3.74 ± 0.07 


log L/L 


0.64 ±0.03 


0.28 ± 0.04 


0.48 ± 0.04 


0.24 ±0.04 


0.41 ±0.03 


0.40 ± 0.03 


BC 


0.01 


-0.03 


0.02 


-0.01 


-0.02 


-0.03 


My 


3.13 ±0.08 


4.08 ± 0.09 


3.52 ±0.10 


4.15 ±0.10 


3.75 ± 0.07 


3.77 ±0.07 


Vo - M v 


6.51 ±0.11 


6.54 ±0.11 


5.83 ±0.13 


5.82 ±0.14 


5.96 ± 0.09 


5.94 ± 0.09 


Distance (pc) 


201 ± 10 


203 ± 10 


147 ±9 


146 ±9 


156 ±7 


154 ±7 


Abundance: 














[Fe/H] 


±0.10 


±0.15 


-0.20 


± 0.12 


-0.27 


±0.07 



t Not corrected for interstellar absorption/reddening. 



peratures, whereas higher interstellar absorption would im- 
ply higher temperatures and therefore a nearly unchanged 
distance. WZ Oph therefore does not formally belong to 
the group of eclipsing binaries within 125 pc discussed by 
Popper ( 1998b ), which could be of use to improve the radia- 
tive flux scale, provided more precise astrometric parallaxes 
were available. 

4.1. Interstellar reddening and absorption 

The interstellar reddenings of AD Boo, VZ Hya, and 
WZ Oph have been determined as follows: 

For AD Boo, we adopt an interstellar reddening of 
E(b — y) = 0.025 ± 0.015, which is derived from the 
calibration by Olsen (|1988[) and the standard uvbyP in- 
dices at phase 0.0, where a significant fraction of the pri- 
mary component is eclipsed. Ideally, indices at phase 0.5, 
where the secondary component is totally eclipsed, should 
be used, but such /? observations are unfortunately not 
available. Combining instead theoretical (3 values (Moon 
& Dworetsky 1 1985ft and the individual uvby indices, red- 
denings of E(b — y) » 0.01 ± 0.02 are obtained for both 
components. Alternatively, 2MASS JHK S photometry for 
stars within 1° of AD Boo on the sky and closer than 200 pc 
show a mean absorption of Ay — 0.18±0.04, corresponding 
to E(b — y) — 0.04 ±0.01 (Knude, private communication). 
Finally, the model by Hakkila et al. (|1997p yields Ay = 0.10 
or E(b — y) = 0.02, in the direction of and at the distance 
of AD Boo, and the Schlegel et al. (1998) dust maps give 
the same amount of total extinction. 

For VZHya, the calibration by Olsen and the stan- 
dard uvbyft indices at phase 0.5, where the secondary is 



nearly totally eclipsed, lead to negligible interstellar red- 
dening, as do the indices at phase 0.0, where a signifi- 
cant fraction of the primary component is eclipsed. On the 
other hand, 2MASS JHK S photometry for stars within 
30' of VZHya on the sky and closer than 150 pc show 
a mean absorption of Ay = 0.17 ± 0.06, corresponding 
to E(b — y) = 0.040 ± 0.014. Furthermore, the model by 
Hakkila et al. also yields Ay = 0.17 in the direction of and 
at the distance of VZ Hya, whereas the Schlegel et al. dust 
maps give a total extinction of Ay — 0.12, corresponding 
to E(b — y) = 0.028. To encompass the different results we 
have adopted E(b -y) = 0.02 ± 0.02. 

For WZ Oph we adopt an interstellar reddening of 
E(b - y) = 0.033 ± 0.012, which is derived from the cal- 
ibration by Olsen and the standard uvbyfi indices out- 
side of eclipse, which represent the nearly identical com- 
ponents well. 2MASS JHK S photometry for stars within 
30' of WZ Oph on the sky and closer than 180 pc show 
a mean absorption of Ay — 0.10 ± 0.04, corresponding 
to E(b - y) = 0.024 ± 0.009. In the direction and at the 
distance of WZ Oph, the model of Hakkila et al. yields a 
very high absorption of Ay = 0.52, whereas the Schlegel 
et al. dust maps give of total extinction of Ay = 0.28 or 
E(b-y) = 0.066. 



5. Discussion 

In the following, we compare the new dimensions obtained 
for AD Boo, VZ Hya, and WZ Oph with properties pre- 
dicted from some of the latest stellar evolutionary models. 
Below, we first briefly present the key ingredients of the 
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selected grids, referring to the original papers for full de- 
scriptions. 

The confrontation between observations and theory for 
these three systems can be constrained better than in most 
previous studies, first of all because our [Fe/H] determina- 
tions and abundance results for other heavy elements are 
probably among the most detailed done to date for binaries. 
This means that we can select models not only with respect 
to mass, but also with respect to Z. On the other hand, for 
a given mass and Z, the observable properties of models 
at a given age such as radius, effective temperature, and 
luminosity still depend on the adopted input physics, in- 
cluding the treatment of core and envelope convection, dif- 
fusion, etc., and on the assumed element mixture and Y,Z 
relation as well. In order to constrain such free model pa- 
rameters, accurate results from many binaries are required, 
supplemented by other results such as those from astero- 
seismology and cluster analysis. We will focus on the three 
binaries studied here; a larger sample will be dealt with in 
forthcoming papers, cf. Fig. [TJ 

From the binary perspective, masses and radii are the 
most direct parameters available, free of any scale de- 
pendent calibrations. The M — R plane is therefore well 
suited for isochrone tests, especially when the binary com- 
ponents are significantly different. In addition, the T e g — R 
or T e ff — log(g) planes allow tests of model temperatures, 
which for a given mass and Z depend, e.g., on the abun- 
dance mixture, Y , core overshoot treatment, and surface 
convection efficiency. In our case, the binary temperatures 
are based on several calibrations and methods, and fur- 
thermore the temperature difference between the compo- 
nents is well constrained from the light curves; cf. Sect. 12.51 
Adopting the usual assumption of coeval formation from 
the same raw material, identical ages must be reached in 
these planes, as well as in the M — log(L) plane. 

5.1. Summary of models 

As a reference, we have adopted the Y 2 (Yonsei-Yale) evo- 
lutionary tracks and isochrones by Demarque et al. (|2T)04lPl. 
which are based on up-to-date input physics and include an 
improved core overshoot treatment where Aos — ^ov/H p 
depends on mass, and which also takes into account the 
composition dependence of M™™ The mixing length 
parameter in convective envelopes is calibrated using the 
Sun, and is held fixed at l/H p = 1.7432. Helium diffusion 
has been included. All tracks are started at the pre-main- 
sequence birth line. Tracks and isochrones are available for 
a broad range of masses, ages, and heavy element abun- 
dances. The enrichment law Y = 0.23 + 2Z is adopted, 
together with the solar mixture by Grevesse et al. ([1996), 
leading to (X,Y,Z) & = (0.71564,0.26624,0.01812). Models 
with enhanced a-element abundances are also included; see 
Kim et al. (|2002p for further details. We have used the abun- 
dance, mass, and age interpolation routines provided by the 
Y 2 group. 

Another recent example of stellar evolution calcu- 
lations is given by the extensive Victoria-Regina grids 



http : / /www . astro . yale . edu/ demarque/yystar . html 

10 Defined as "the mass above which stars continue to have 
a substantial convective core even after the end of the pre-MS 
phase." 



(VandenBerg et al.,[5(]O60, which adopt a somewhat dif- 
ferent physically based description of the core overshoot 
with F over depending on mass and metallicity and cali- 
brated observationally via cluster CMD's. Diffusive pro- 
cesses are not treated. All tracks have been computed from 
the pre-main-sequence stage. The VRSS grids (scaled-solar) 
used here adopt Y = 0.23544 + 2.2Z and the solar mix- 
ture by Grevesse & Noels flUM]), leading to (X,Y,Z) Q 
= (0.7044,0.2768,0.0188). The mixing length parameter in 
convective envelopes is calibrated using the Sun and held 
at otMLT = 1-90. For further details on the chemical abun- 
dances adopted for the VR0A, VR2A, and VR4A grids 
(different [a/Fe]) we refer the reader to VandenBerg et al. 
(2000}. We have used the isochrone interpolation routines 
provided by the Victoria-Regina group; track interpolations 
have not been done. 

A third source of modern evolutionary tracks and 
isochrones for specific masses and ages is the BaSTI 
database (Pietrinferni et al., I2004[P^1 , available for sev- 
eral compositions. It is based on a recent version of the 
FRANEC evolutionary code (Degl'Innocenti et al.. |2007|) . 
Models without and with core overshoot are available; in 
the latter case the choice of X ov depends on mass but not 
on metallicity. Y = 0.245 + 1.4Z and the solar mixture by 
Grevesse & Noels (|1993|) are used, and the mixing length 
parameter in convective envelopes is calibrated using the 
Sun and held fixed at l/H p — 1.913. Mass loss is included, 
but atomic diffusion is taken into account only in the stan- 
dard solar model. All tracks have been computed from the 
pre-main-sequence phase. 

A comparison of the main sequence models between 1.1 
and 1.5 Mq for identical [Fe/H] values close to those of our 
stars reveals some differences between the three grids de- 
scribed above. This is exemplified in Figs. [TBI and [P9l In the 
lower main-sequence band, the Y 2 and the Victoria-Regina 
tracks agree quite well. The isochrones are slightly shifted, 
and for part of the mass interval the effect increases in the 
upper part of the main-sequence band due to differences in 
the treatment of core overshoot and the ramping procedure 
adopted. The BaSTI tracks differ somewhat in shape and 
location, and consequently significant differences arc seen 
in the isochrones as well. 

Before focusing on the confrontation between binary ob- 
servations and models, we point out that we have not con- 
sidered here models based on scaling of the most recent 
solar mixture (Grevesse et al. I2007|) . which implies a sig- 
nificantly lower metallicity of Zq w 0.013 and a different 
conversion between [Fe/H] and Z. At present, solar models 
computed with these abundances disagree with precise re- 
sults inferred from helioseismology, and the cause(s) have 
yet to be identified (e.g. Basu & Antia l20"0"7j) . 

5.2. Comparison with the binary results 

For AD Boo, comparisons between the absolute dimensions 
and Y 2 models calculated for the measured [Fe/H] are 
shown in Figs. [2014221 As seen in Fig. [211 an isochrone cor- 
responding to ~1.75 Gyr would fit both components very 
well, and for the primary component this age is also de- 
duced from Figs. l20l and l22l The abundance uncertainty of 



|http : //wwwl . cadc-ccda. hia-iha.nrc-cnrc . gc . ca/cvo, 
community/VictoriaReginaModels/ 
^ |http : //www . te . astro . it/BASTI/index.php| 
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Fig. 18. Comparison between Y 2 (red) and Victoria- 
Regina VRSS (blue) models for [Fe/H] = -0.190 and 
[a/Fe] = 0.0. 

Y 2 : (XXX) = (0.73400,0.25400,0.01200). VRSS: (XX, Z) 
= (0.72456,0.26294,0.01250). Tracks (solid lines) for 1.1, 
1.2, 1,3. 1.4, and 1.5 Mq and isochrones (dashed) between 
0.5 and 2.5 Gyr (step 0.5 Gyr) are shown. 
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Fig. 19. Comparison between Y 2 (red) and BaSTI core 
overshoot (blue) models for [Fe/H] = —0.253 and [a/Fe] = 
0.0. 

Y 2 : (XXX) = (0.73865,0.25090,0.01045). BaSTI: 
(XXX) = (0.7310,0.2590,0.0100). Tracks (solid lines) for 
1.1, 1.2, 1.3, 1.4, and 1.5 Mq and isochrones (dashed) 
between 0.5 and 2.5 Gyr (step 0.5 Gyr) are shown. 



±0.15 dex translates to an age uncertainty of about ±0.15 
Gyr, comparable to that from mass and radius uncertain- 
ties. The model tracks fit the components well but predict a 
~ 100 K lower temperature difference between the compo- 
nents than established from both photometry and the em- 
pirical flux scale; see Sect. 0J Consequently, the predicted 



Fig. 20. AD Boo compared with Y 2 models for (XXX) 
= (0.7028,0.2748,0.0224), equivalent to the measured abun- 
dance [Fe/H]= ±0.10 ±0.15 for [a/Fe] = 0.0. Tracks for 
the component masses (solid lines) and isochrones from 0.5 
to 3.0 Gyr (dashed; step 0.5 Gyr) are shown. The uncer- 
tainty in the location of the tracks coming from the mass 
errors are indicated (dotted lines). To illustrate the effect 
of the abundance uncertainty, tracks for [Fe/H] = +0.25 
and [a/Fe] = 0.0 (dash-dot), corresponding to (XXX) = 
(0.6785,0.2910,0.0305), are included. 

model luminosity of the secondary is also marginally higher 
than observed (Fig. |2"2")) . If real, the slight discrepancy may 
be related to the activity of the secondary component. The 
same picture is seen for Victoria-Regina models, which give 
a slightly lower age of ~1.5 Gyr. On the other hand, the 
BaSTI models fit AD Boo less well. For the observed masses 
and radii, they predict quite different ages of about 2.25 
and 1.75 Gyr for the primary and secondary, respectively. 
This is probably related to the different shape of the BaSTI 
tracks (Fig. [19]) in the lower main-sequence region. 

Figs. l2"3Tf2l)l compare the absolute dimensions for 
VZ Hya with the Y 2 models calculated for the measured 
[Fe/H]. As seen in Fig. [24] an isochrone for an age of 1.25 
Gyr fits both components within errors, although a slightly 
steeper isochrone slope would formally improve the fit still 
more. The abundance uncertainty of ±0.12 dex corresponds 
to about ±0.25 Gyr in age, similar to the contribution from 
mass and radius uncertainties. The model tracks are slightly 
hotter than observed, and consequently the ages derived 
from Figs. |2"31 andl2"5"lbecome too high and low, respectively. 
Both tracks would fit well for a slightly higher [Fe/H] of 
—0.15, safely within the uncertainty of the measured value. 
Adopting [Fe/H] = —0.15, identical ages of 1.25 Gyr are ob- 
tained from all three diagrams. Victoria-Regina tracks for 
the observed [Fe/H] fit the components well, being slightly 
cooler than the Y tracks. For the observed masses and 
radii, ages become lower, about 1.0 and 0.7 Gyr for pri- 
mary and secondary, respectively; they agree within errors. 
The BaSTI models yield about 1.25 and 0.75 Gyr, and in 
this case the difference is significant. 

Whereas the Y 2 models in general match the observed 
properties of AD Boo and VZHya well, the situation for 
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Fig. 21. AD Boo compared with Y 2 models for (X,Y,Z) 
= (0.7028,0.2748,0.0224), equivalent to [Fe/H] = +0.10 for 
[a/Fe] = 0.0. Isochrones from 0.5 to 3.0 Gyr (step 0.5 
Gyr) are shown. The 1.5 Gyr isochrone for [Fe/H] = +0.25 
and [a/Fe] = 0.0 (dashed), corresponding to (X,Y,Z) = 
(0.6785,0.2910,0.0305), is included for comparison. 




Fig. 22. AD Boo compared with Y 2 models for (X,Y,Z) 
= (0.7028,0.2748,0.0224), equivalent to [Fe/H] = +0.10 for 
[a/Fe] = 0.0. Isochrones from 0.5 to 3.0 Gyr (step 0.5 Gyr) 
are shown. 

WZ Oph is completely different, as seen from Fig. [26] 
Models for the measured [Fe/H] are much too hot, both 
for [a/Fe]= 0.0 and +0.1 (cf. Sect. EHH)) . and the same 
difference is seen for Victoria-Regina and BaSTI models. 
In order to reproduce WZ Oph for the observed [Fe/H], an 
unrealistically high [a/Fe] of about +0.5 would be required. 
Alternatively, the models would need to have a significantly 
lower helium content and/or an envelope mixing length pa- 
rameter well below l/H p — 1.7432. The latter can not a 
priori be ruled out, since the photometry and spectroscopy 
indicate the presence of some degree of surface activity 



Fig. 23. VZ Hya compared with Y 2 models for (X,Y,Z) 
= (0.73478,0.25348,0.01174) equivalent to the measured 
abundance [Fe/H] = -0.20 ± 0.12 for [a/Fe] = 0.0. Tracks 
for the component masses (solid lines) and isochrones from 
0.5 to 2.5 Gyr (dashed; step 0.5 Gyr) are shown. The uncer- 
tainty in the location of the tracks coming from the mass 
errors are indicated (dotted lines). To illustrate the effect 
of the abundance uncertainty, tracks for [Fe/H] = —0.08 
and [a/Fe] = 0.0 (dash-dot), corresponding to (X,Y,Z) = 
(0.72425,0.26050,0.01525), are included. 



on both components of WZ Oph; see Sect. 13.3.21 It would 
therefore be of great interest to have available models for 
the observed [Fe/H] and a range of l/H p values. Such mo- 
dels have been published by Baraffe et al. (|1998P^1 but only 
for solar metallicity. They suggest that l/H p ~ 1 or even 
lower would be needed, although models computed speci- 
fically for this case would be required in order to evaluate 
the combined effects on the stellar radii and temperatures 
in detail. The possible connection between activity and con- 
vection, as indicated from studies of solar-type binaries, has 
been discussed by Clausen et al. ( 199jl|), Torres et al. ([2006), 
and others. 

As stated in Sect. 14.11 we have adopted an interstellar 
reddening of E(b - y) = 0.033 ± 0.012 for WZ Oph, equi- 
valent to Ay — 0.14 ± 0.05 for a standard extinction law. 
There were, however, some indications of higher absorption 
from other sources. If we therefore somewhat arbitrarily 
double it and assume Ay = 0.28 and E(b — y) = 0.066, 
the effective temperatures derived from the uvby indices 
increase by about 250 K, i.e. to about 6400 K. This has 
an effect on the abundances. From the Fe I lines we would 
derive [Fe/H] ~ —0.13 for both components, whereas the 
results from the Fe II lines would be practically unchanged 
at —0.27. The photometric determinations from the cali- 
bration by Holmberg et al. (2007) would increase by about 
0.12 dex for both components. The discrepancy between 
neutral and ionized Fe lines argues against the higher tem- 
peratures, but assuming them, together with the higher 
abundance from the Fe I lines, does on the other hand re- 
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Fig. 24. VZ Hya compared with Y 2 models for (X,Y,Z) 
= (0.73478,0.25348,0.01174) equivalent to [Fe/H] = -0.20. 
Isochrones from 0.5 to 2.5 Gyr (step 0.5 Gyr) arc 
shown. The 1.5 Gyr isochrone for [Fe/H] = —0.08 and 
[a/Fe]= 0.0 (dashed), corresponding to (X,Y,Z) = 
(0.72425,0.26050,0.01525), is included for comparison. 
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Fig. 25. VZ Hya compared with Y 2 models for (X,Y,Z) 
= (0.73478,0.25348,0.01174) equivalent to [Fe/H] = -0.20. 
Isochrones from 0.5 to 2.5 Gyr (step 0.5 Gyr) are shown. 



move the model discrepancies, as seen in Fig. [57J The age 
of WZ Oph then becomes about 3 Gyr. 

An independent check would be to derive intrinsic 
colours and/or effective temperatures from line-depth ra- 
tios measured in the disentangled FEROS spectra, see e.g. 
Gray (|1994p . but WZOph is unfortunately slightly too hot 
for the calibrations available. Adopting a subset of those 
by Kovtyukh et al. (j2"Ut?3"|) . valid to 6100 or 6150 K, we de- 
rive mean effective temperatures close to 6200 K for both 
components. The individual results from the 9 (primary) 
and 16 (secondary) measured line-depth ratios scatter by 
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Fig. 26. WZ Oph compared with Y 2 models for (X,Y,Z) 
= (0.73985,0.25010,0.01005), equivalent to the measured 
abundance [Fe/H] -0.27 ±0.07 for [a/Fe] = 0.0. Tracks for 
the component masses (solid lines) and isochrones from 0.5 
to 3.0 Gyr (dashed; step 0.5 Gyr) are shown. The uncer- 
tainty in the location of the primary track coming from the 
mass error is indicated (dotted line). For comparison, tracks 
for (X,Y,Z) = (0.7349,0.2534,0.0117) and [a/Fe] = +0.1, 
also equivalent to [Fc/H] = —0.27, are included (dash-dot). 




6800 6600 6400 6200 6000 
T e „ [K] 



Fig. 27. Assuming E(b-y) = 0.066: WZ Oph compared 
with Y 2 models for (X,Y,Z) = (0.72899,0.25734,0.013670), 
equivalent to [Fe/H] = -0.13 for [a/Fe] = 0.0. Tracks for 
the component masses (full drawn) and isochrones from 
0.5 to 3.5 Gyr (dashed; step 0.5 Gyr) are shown. See text 
for details. 

100-200 K, and since we are at the edge of the calibrations, 
the method does not at the moment provide new firm in- 
formation. 

Although we cannot pinpoint the reason for the discrep- 
ancy between the observations and the models for WZ Oph, 
we suspect they are connected to composition (He) or con- 
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vection anomalies in the stars and/or the determination of 
interstellar reddening, rather than to physical or numerical 
problems with the model calculations. This is supported by 
the general good agreement seen for AD Boo and especially 
VZ Hya, but detailed studies of additional systems similar 
to these are needed for full clarification. Work on several 
new F-G type binaries is in progress. 

6. Summary and conclusions 

We have presented precise absolute dimensions and abun- 
dances for the three F-type double-lined eclipsing bina- 
ries AD Boo, VZ Hya, and WZ Oph. The results, which are 
based on state-of-the-art analyses of new photometric and 
spectroscopic observations, were used to test the recent Y 2 , 
Victoria- Regina, and BaSTI stellar evolutionary models. 

Masses and radii precise to 0.5-0.7% and 0.4-0.9% have 
been established for the binary components. This level can 
only be reached if excellent observations are available and 
are analysed carefully using adequate tools. Special care 
was taken to identify and avoid possible systematic errors. 
Several methods were applied in order to determine reli- 
able and accurate (100-150 K) effective temperatures and 
interstellar extinctions, still a challenging task whose im- 
portance can not be overestimated. 

System [Fe/H] abundances were based on numerous 
lines and were derived relative to the Sun. We have reached 
precisions between 0.07 and 0.15 dex. Even having many 
Fe I and Fe II lines of different strength and excitation po- 
tential available, the strong correlation between abundance 
and effective temperature is difficult to break. T e ff and mi- 
croturbulence uncertainties contribute significantly to the 
final abundance precision. Abundances for other heavy el- 
ements, based on fewer lines, were also derived; in one case 
(WZ Oph) indications of a slight a-element overabundance 
is seen. Comparing the results for the three binaries clearly 
demonstrates that analyses of disentangled spectra are far 
superior to those of composite spectra, essentially because 
many more lines become available and blending becomes 
less of an issue. 

AD Boo has components with the greatest difference 
in mass and radius among well-studied F-type binaries 
(1.41 Mq, 1.61 R Q ; 1.21 Mq, 1.22 Rq) and an abundance 
of [Fc/H] = +0.10. Slight surface activity is present on the 
secondary component. The system is fitted well by Y 2 and 
Victoria-Regina evolutionary models for identical compo- 
nent ages (as determined from masses and radii) of 1.75 
and 1.5 Gyr, respectively. The fact that the secondary ap- 
pears marginally cooler than the corresponding evolution- 
ary track may be related to mild surface activity. BaSTI 
models predict different ages of 2.25 and 1.75 Gyr for the 
components, and therefore do not fit AD Boo as well. 

VZHya has somewhat different components (1.27 Mq, 
1.31 R Q ; 1.14 Mq, 1.11 Rq) and an abundance of [Fe/H] = 
—0.20. No signs of surface activity are seen. Y 2 and 
Victoria-Regina evolutionary tracks fit VZ Hya well. Nearly 
identical component ages (as determined from masses and 
radii) of 1.25 Gyr are obtained from the Y 2 isochrones, 
whereas the Victoria-Regina calculations give 1.0 and 0.7 
Gyr for the primary and secondary. An even larger differ- 
ence, 1.25 versus 0.75 Gyr, is found when using the BaSTI 
models. 

WZ Oph consists of two nearly identical components 
(1.22 Mq, 1.41 Rq). The [Fe/H] abundance is -0.27, and 



as mentioned above a slight a-element overabundance is 
possible. Photometry and spectroscopy suggest low-level 
surface activity on both components. Contrary to the case 
for AD Boo and VZ Hya, neither of the three model grids are 
able to fit the observations for WZ Oph; models are signif- 
icantly hotter than observed. We tentatively conclude that 
this is caused either by WZ Oph anomalies such as a low 
He content or decreased envelope convection, or by under- 
estimated interstellar reddening. Further observations will 
be required to resolve the issue. 

We believe our study is among the most detailed car- 
ried out to date for main-sequence eclipsing binaries, and 
that it sets new standards for critical tests of stellar evolu- 
tionary models based on accurate binary data. With only 
three systems available, it is premature to draw firm con- 
clusions and suggest specific model shortcomings and/or 
improvements. We find, however, better agreement with 
the Y 2 and Victoria-Regina models than with the BaSTI 
models. In contradiction to this, Tomasella et al. (2007|) re- 
cently found that dedicated BaSTI models, including core 
overshoot and both helium and heavy element diffusion, 
represent the components of the young ~ 1.25Mq system 
V505 Per well. Their study is based on new absolute dimen- 
sions of extraordinary high formal precision. 

As shown in Fig. [H several additional F-type binaries 
are available, and we are presently working on re-analyses 
for some of them, including abundance determinations. We 
are also conducting a larger program on new F-G type sys- 
tems exhibiting various levels of surface activity (Clausen 
et al. I200ip . and we expect some very useful insights to 
come from these studies, and from a parallel program on 
eclipsing binaries in clusters spanning a wide range in age 
and metallicity (e.g. Grundahl et al. I2006p . 
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Appendix A: Radial velocity tables 
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Table A.l. Radial velocities (corrections applied) of AD Boo and resid- 
ual from the final spectroscopic orbit. 



HJD 


Phase 


RV P 


RV S 


(O - C 


)p 


in - c) 

\^ ^ )s 


-2400 000 






kni s 


-l 


km s 


-l 




■ l 




47555. 


.8293 





.3938 


-83. 


,75 


58. 


.17 


-1 


.57 


-2.09 


47598. 


.8459 





.1868 


-114 


.44 


97. 


.36 


-0 


.05 


-0.57 


47612. 


6969 





.8820 


56. 


,10 


-103 


.17 


0. 


,92 


-2.81 


47613. 


.8580 





.4432 


-53 


,94 


25. 


.71 


-0 


.34 


-1.13 


47628. 


.7443 





.6388 


66. 


,83 


-113 


.43 


2, 


.09 


-1.89 


47640. 


.6739 





.4052 


-77. 


,17 


57. 


.34 


-1 


.11 


4.23 


47642. 


.7577 





.4125 


-69. 


.91 


49 


.99 


2, 


.09 


1.63 


47643. 


.6147 





.8267 


75. 


.39 


-124 


.58 


-2 


.14 


1.91 


47644. 


.7409 





.3711 


-93 


.70 


76. 


.77 


-0 


.31 


3.39 


47675. 


.5822 





.2788 


-123 


.92 


103 


.72 


-3 


.01 


-1.84 


47676. 


.5584 





.7507 


88. 


.38 


-142 


.11 


-1 


.23 


-1.48 


47693. 


.7772 





.0738 


-66 


.80 


36 


,38 


-2 


.83 


-2.59 


47700. 


.6648 





.4030 


-77. 


.85 


59 


,20 


-0 


.59 


4.69 


47701. 


.7853 





.9446 


20. 


,07 


-53 


.80 


0, 


,42 


5.02 


47702. 


.6453 





.3603 


-97. 


,54 


79. 


.78 


0. 


,61 


0.84 


47702. 


.7425 





.4073 


-74. 


.63 


52. 


.97 


0. 


,25 


1.24 


47703. 


.6426 





.8424 


70. 


.23 


-120 


.29 


-1 


.99 


0.00 


47703. 


.7439 





.8914 


51. 


,04 


-97. 


.68 


0. 


,62 


-2.88 


47729. 


.5844 





.3819 


-86. 


,88 


72 


,30 


1, 


,36 


4.95 


47731. 


.5528 





.3334 


-106 


.61 


91. 


,39 


1, 


,81 


0.44 


47734. 


.5463 





.7803 


88 


,81 


-141 


.47 


1, 


,12 


-3.09 


47957. 


.8075 





.6982 


81. 


,97 


-130 


.23 


-2 


.07 


3.88 


47958. 


.8221 





.1886 


-115 


.20 


94. 


.94 


-0 


.34 


-3.53 


47959. 


.8666 





.6935 


81. 


,81 


-127 


.99 


-1 


.19 


4.90 


47960. 


.9249 





.2050 


-119 


,65 


105 


.59 


-1 


.20 


2.92 


48279. 


.8560 





.3669 


-96. 


11 


72. 


.02 


-0 


.82 


-3.58 


48280. 


.9257 





.8840 


53. 


,24 


-104 


.72 


-0 


.95 


-5.51 


48311 


.9322 





.8716 


59. 


,11 


-101 


.51 


-1 


.02 


4.64 


48499. 


.5443 





.5577 


19. 


,69 


-53 


,87 


-1 


.43 


6.67 


48638. 


.8910 





.9138 


36 


,51 


-79 


,95 


-1 


.70 


0.57 


48639. 


.9724 





.4365 


-57. 


,70 


32 


,06 


0. 


.06 


0.35 


48699. 


.9584 





.4319 


-61. 


.66 


32 


,15 


-1 


.12 


-2.80 


48701. 


.8654 





.3537 


-101 


.73 


83. 


.27 


-0 


.83 


1.12 


48704. 


,9191 





.8298 


79. 


,57 


-125 


.73 


3, 


,02 


-0.37 


48705. 


.8840 





.2962 


-120 


.74 


107 


.80 


-2 


.53 


5.40 


48706. 


.8287 





.7528 


92. 


,19 


-139 


.10 


2, 


,59 


1.51 


48711. 


.8041 


o 


.1578 


-106 


.84 


83 


,82 


-1 


.50 


-3.52 


48725. 


.8331 





.9390 


21. 


,13 


-60. 


,51 


-2 


.04 


2.42 


48728. 


.8111 





.3785 


-93 


,65 


70. 


,35 


-3 


.74 


1.05 


48730. 


.9634 





.4188 


-69. 


,17 


44. 


,35 


-0 


.84 


0.28 


48754. 


.5676 





.8284 


77. 


,25 


-127 


.14 


0. 


.26 


-1.27 


48759. 


.7935 





.3545 


-99 


,14 


87. 


.17 


1, 


,47 


5.36 


48783. 


.6614 





.8915 


52. 


,54 


-97 


,21 


2. 


.17 


-2.48 


48783. 


.7736 





.9457 


18. 


,44 


-58. 


.86 


-0 


.54 


-0.83 


48784. 


.7101 





.3984 


-78. 


,44 


62. 


.33 


1, 


,33 


4.89 


48810. 


.6101 





.9177 


34. 


,75 


-77. 


.67 


-1 


.21 


0.21 


48824. 


.6164 





.6879 


79. 


,52 


-130 


.74 


-2 


.12 


0.57 


48828. 


.5744 





.6011 


44. 


.50 


-88, 


.00 


-1 


.96 


2.16 


48840. 


.5419 


o 


.3858 


-87. 


,71 


69 


.26 


-1 


.42 


4.19 


49106. 


.8924 





.1318 


-91. 


,82 


76, 


.57 


2, 


,87 


1.68 


49111. 


.6651 





.4388 


-57. 


,14 


30 


.24 


-0 


,78 


0.17 


49112 


.6802 





.9294 


27. 


39 


-73 


.06 


-1 


.63 


-3.29 


49114 


.7754 





.9422 


22. 


.66 


-59, 


.87 


1, 


.46 


0.75 


49115 


.7720 





.4239 


-65. 


.94 


41, 


.63 


-0 


.59 


1.05 


49117 


,7719 





.3906 


-84. 


.69 


61 


.93 


-0 


.83 


-0.30 


49136. 


.7641 





.5709 


28. 


.22 


-72. 


.23 


-0 


.97 


-2.26 


49138. 


.8295 





.5692 


24. 


,31 


-67. 


.62 


-3 


.89 


1.19 


49140. 


.6382 





.4435 


-54. 


.56 


19. 


.20 


-1 


.14 


-7.43 



J.V. Clausen et al.: AD Boo, VZHya, and WZOph 
Table A.l. continued. 



HJD Phase RV P RV S (O - C) p (O - C) 



-2400 000 






km s 


-1 


km s 


-1 


km s 


■ 1 


km s 


-1 


4yi41. 


a c a 1 
.5541 





.9345 


28, 


.08 


-65 


.19 


2. 


,17 


0. 


.94 


A H1 A o 


bozo 





.4315 


-59. 


.39 


41. 


.82 


1, 


,38 


6. 


.59 


/i n 1 /ic 

4yi4o 


7A70 

.7y7z 





,420b 


-66. 


,61 


40. 


,58 


0. 


,71 


-2 


.30 


A n 1 70 

4yi To. 


.olio 





.3820 


-85. 


.69 


65. 


,83 


2. 


,53 


-1 


.49 


4yi74. 


.5980 





.8587 


65. 


,86 


-114 


.47 


0. 


,03 


-1 


.0b 


4yi T4. 


75yu 





.9355 


25. 


,50 


-66. 


.53 


0. 


,77 


-1 


.78 


4yi7o. 


bbUy 





.3724 


-93 


.04 


71. 


.89 


-0 


.27 


-0 


,7b 


4y4(Jl 







.5718 


76. 


,53 


-127 


.27 


-0 


.54 


-1 


.31 


AC\A AT 

4y44 f . 


7AA7 

. / 00 / 





,8b84 


63 


,10 


-116 


.37 


1, 


.52 


-8. 


.53 


4y5D4. 


Ct ACi 

.b!4y 





.3813 


-86. 


.88 


69. 


.00 


1, 


,68 


1, 


.28 


A C\TCl 1 


0070 

.oo7z 





.7127 


85. 


,16 


-131 


.47 


-1 


.55 


5, 


.70 


4y mo. 


ooOl 





.b524 


69 


,42 


-116 


.29 


-0 


,85 


1, 


.72 


4y7t)7. 


.ouyu 





.5993 


45. 


,62 


-89. 


.81 


0. 


,15 


-0 


.81 


/I noo c 


071 O 

.o71z 





.bb48 


74. 


,12 


-122 


.47 


-0 


,65 


0, 


.80 


4yo4o. 


001 a 

.©zlb 





.3415 


-104 


,74 


86 


.98 


0. 


.85 


-0 


.Ob 


/i no to 
4yo7o. 


T A A O 

. /44z 





.8052 


83. 


.09 


-125 


.64 


-0 


.21 


7, 


.00 


4ysyi. 


.bzll 





.4454 


-52. 


,53 


26 


,95 


-0 


.92 


2, 


.43 


49901. 


.5946 





.2b73 


-123 


.33 


105 


,21 


-1 


.30 


-1 


.05 


a nnnn 

4yyuy. 


c no 1 





.1330 


-96. 


,19 


77. 


,05 


-0 


.94 


1 


.51 


4yyz(j. 







.4487 


-48. 


.66 


24. 


,62 


1. 


,49 


1 


.81 


a nnon 

49939. 


C 7AO 

570z 





.b235 


60. 


.30 


-99 


,67 


2. 


,46 


3. 


.80 


49955 


ri /in 





.1544 


-105 


.60 


96. 


,18 


2. 


,08 


6. 


.10 


49979 


cm o 

.501o 





.9250 


36 


,07 


-70. 


,40 


4, 


,43 


2. 


.43 


50119. 


.obby 





.7736 


87. 


,83 


-134 


.44 


-0 


,62 


4, 


.82 


50126. 


m a n 





.1800 


-113 


.41 


96. 


.07 


-0 


.86 


0. 


.30 


50140. 


7AAA 

7yzy 





.8886 


51. 


.90 


-99 


.18 


0, 


.06 


-2 


.73 


5014b 


.OOZZ 





.8320 


77. 


.74 


-125 


.27 


1, 


,91 


-0 


.76 


50154. 


.oolb 





.6745 


76. 


.19 


-127 


.41 


-1 


.71 


-0 


.48 


501 ( 







.4116 


-73 


.04 


54. 


.00 


-0 


.57 


5. 


.10 


5050 r 


O A 7H 

.84/0 





.3117 


-113 


.80 


96. 


.01 


0, 


.98 


-2 


.37 


50514. 


001 7 





.6879 


80. 


,94 


-130 


.93 


-0 


.70 


0. 


.37 


5051b. 


AA/J 1 

yubi 





.6906 


85. 


,91 


-131 


.11 


3 


.60 


0, 


.97 


50534. 


OO/IO 

.504z 





.3565 


-100 


,21 


76 


,31 


-0 


.45 


-4 


.51 


50535 


.b4L)y 





.7464 


89. 


,82 


-139 


.25 


0. 


.23 


1, 


.35 


50535 


.bobl) 





.7537 


91. 


,99 


-146 


,77 


2. 


.40 


-6 


.18 


50535. 


b710 





.7610 


90 


,27 


-135 


.81 


0. 


,91 


4, 


,52 


C O C /I o 

50540. 


.ozOb 





.2502 


-124 


.08 


109 


.00 


-1 


.43 


1, 


,41 


50541. 


7Aon 

jy^y 





.7206 


86. 


,82 


-140 


.28 


-0 


.99 


-1 


.75 


50543 


.oobo 





.7080 


83. 


.48 


-131 


.91 


-2 


.45 


4, 


,41 


50561 


7Aon 





.3469 


-105 


.72 


80 


,35 


-2 


.15 


-4 


.92 


50563 


71 i r 

.7115 





.3149 


-113 


.48 


96. 


,44 


0. 


,46 


-0 


.9b 


50554. 


/JAAA 

byzy 





.7893 


86. 


,35 


-139 


.03 


-0 


.04 


-2 


.17 


50555. 


b50o 





.2523 


-120 


.05 


101. 


,11 


2. 


,59 


-0 


.47 


50567 


171 A 

.7134 





.2493 


-123 


,07 


104. 


.91 


-0 


.42 


-2 


.68 


50568. 


.by5y 





.7243 


90 


,61 


-138 


.11 


2. 


,38 


0. 


.90 


50559 


TTT(\ 

.777y 





.2473 


-121 


,60 


109 


.17 


1, 


,04 


1, 


,60 


50571 


A 1 

.5541 





.1590 


-105 


.86 


88 


,54 


-0 


,09 


0. 


.69 


50573. 


Toon 
. /839 





.1836 


-116 


.28 


98. 


,59 


-2 


.72 


1 


.63 


50589. 


.5959 





.8267 


77. 


,20 


-127 


.53 


-0 


.33 


-1 


.03 


50590. 


5390 





.3309 


-107 


.69 


93 


.36 


1, 


,55 


1 


.40 


50591. 


.bb85 





.8285 


75. 


,92 


-124 


.18 


-1 


.04 


1 


.05 


50592. 


.6487 





.3023 


-116 


,76 


100. 


,03 


0. 


.21 


-0 


.91 


50593. 


.5959 





.7602 


87. 


,43 


-140 


.95 


-1 


.97 


-0 


.57 


50595. 


.5517 





.7587 


86. 


,84 


-135 


.87 


-2 


.62 


4. 


.57 


5059b. 


.5967 





.2107 


-117 


,38 


102. 


.96 


2. 


,05 


-0 


.80 


50597. 


.b074 





.6992 


82. 


,68 


-135 


.60 


-1 


.58 


-1 


.24 


50519. 


,b0b4 





.3329 


-106 


.61 


94 


,13 


1, 


,98 


2. 


.99 


50519. 


.b215 





.3402 


-106 


.01 


91 


,36 


0. 


.06 


3. 


.16 


50520. 


.6194 





.8225 


79. 


.11 


-126 


.68 


0. 


.33 


1. 


.28 



J.V. Clausen et al.: AD Boo, VZHya, and WZOph 
Table A.l. continued. 



11.) 1 J 


Phase 


pi/ 

HVp 


KV S 


(n n\ 


\y - ^)s 


— 1 4UU UUU 




kms -1 


kms -1 


kms -1 


km s _1 


50621.5810 


0.2873 


-120.35 


107.51 


-0.60 


3.32 


50622.6012 


0.7805 


89.09 


-133.20 


1.41 


5.16 


50623.6005 


0.2635 


-122.58 


108.27 


-0.31 


1.13 


50624.6042 


0.7487 


87.14 


-139.68 


-2.47 


0.94 


50625.6138 


0.2367 


-122.61 


108.53 


-0.33 


1.37 


50627.6490 


0.2204 


-123.42 


109.14 


-2.59 


3.68 


50629.5948 


0.1610 


-106.47 


90.60 


0.01 


1.92 
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Table A. 2. Radial velocities (corrections applied) of 
VZHya and residuals from the final spectroscopic orbit. 



H.7D 


Phase 


RV P 


RV 3 


(O - C% 


(O - C). 


-2 400 000 






km s~ 1 


km s~ 1 


km s 1 


km s~ 1 


47555 


.7226 





.8107 


82.13 


-101.37 


-1.40 


0.94 


47558 


.7158 





.8413 


75.81 


-89.61 


0.66 


3.40 


47568 


.6242 





.2530 


-99.07 


101.23 


0.40 


0.51 


47570 


.7110 





.9715 


13.91 


-20.42 


1.56 


2.92 


47574 


.6105 





.3141 


-93.26 


92.10 


-1.38 


-0.20 


47575 


.6910 





.6862 


82.32 


-100.78 


-0.50 


0.74 


47577 


.6268 





.3527 


-79.70 


81.21 


0.70 


1.64 


47583 


.6110 





.4132 


-54.33 


49.75 


-0.51 


-0.33 


47587 


.6242 





.7950 


86.63 


-105.73 


0.05 


-0.03 


47603 


.5188 





.2678 


-98.83 


100.09 


0.06 


0.00 


47609 


.5595 





.3477 


-80.89 


83.10 


1.27 


1.58 


47613 


.5894 





.7352 


90.37 


-110.66 


0.43 


-1.24 


47628 


.5692 





.8930 


54.76 


-71.57 


0.23 


-1.44 


47642 


.5540 





.7082 


85.75 


-106.46 


-1.35 


-0.19 


47645 


.5312 





.7333 


89.68 


-108.31 


-0.15 


0.99 


47664 


.5534 





.2830 


-97.70 


100.29 


-0.25 


1.81 


47670 


.6171 





.3708 


-74.23 


70.09 


-0.82 


-1.72 


47671 


.6217 





.7167 


87.91 


-106.94 


-0.37 


0.64 


47672 


.6262 





.0626 


-39.66 


38.82 


1.29 


3.02 


47899 


.7202 





.2549 


-101.14 


100.79 


-1.70 


0.10 


47905 


.7028 





.3149 


-92.07 


93.39 


-0.36 


1.27 


47928 


.7245 





.2416 


-98.52 


102.25 


0.83 


1.66 


47941 


.6387 





.6882 


83.84 


-100.78 


0.56 


1.25 


47956 


.6250 





.8482 


72.18 


-91.00 


-0.66 


-0.55 


47957 


.6126 





.1883 


-91.30 


96.46 


1.14 


3.54 


47959 


.5793 





.8654 


65.89 


-83.00 


-0.56 


0.36 


47960 


.5813 





.2105 


-95.81 


98.29 


0.76 


0.78 


48191 


.8651 





.8454 


74.21 


-89.36 


0.41 


2.16 


48191 


.8800 





.8505 


72.28 


-87.06 


0.25 


2.50 


48226 


.7915 





.8712 


64.13 


-83.69 


-0.02 


-2.88 


48230 


.9468 





.3019 


-96.99 


89.83 


-2.51 


-5.36 


48236 


.8621 





.3386 


-87.41 


83.89 


-2.27 


-0.94 


48252 


.8201 





.8332 


77.88 


-94.81 


0.22 


0.99 


48256 


.8706 





.2279 


-98.67 


99.73 


-0.10 


0.00 


48258 


.8791 





.9195 


42.00 


-56.01 


0.56 


-0.40 


48259 


.7319 





.2131 


-96.77 


97.61 


0.17 


-0.31 


48264 


.7826 





.9521 


22.39 


-34.46 


-1.16 


1.30 


48280 


.7728 





.4578 


-30.24 


23.96 


-0.82 


0.95 


48281 


.7752 





.8030 


82.86 


-103.45 


-2.28 


0.64 


48290 


.6911 





.8729 


63.96 


-84.32 


0.52 


-4.30 


48315 


.6410 





.4636 


-24.29 


18.58 


1.82 


-0.76 


48354 


.5562 





.8627 


68.11 


-84.77 


0.59 


-0.22 



Table A. 3. Radial velocities (corrections applied) of 
WZOph and residuals from the final spectroscopic orbit. 



H.ID 


Phase 


RV P 


RV S 


(O - C) p 


(0-C*) s 


-2 400 000 






km s — 1 


km s - 1 


km s 


1 


km s 1 


47583. 


9253 


0. 


4058 


-78.81 


22.61 


-0 


.94 


1.19 


47605. 


8830 


0. 


6545 


45.26 


-101.95 





.37 


0.06 


47608. 


8786 





3705 


-94.35 


35.97 


-1. 


.48 


-0.54 


47612. 


.8755 





3259 


-107.37 


49.40 


-0 


.15 


-1.54 


47628. 


8294 


0. 


1394 


-95.42 


41.35 


1. 


.13 


1.14 


47639. 


8265 





.7681 


58.80 


-117.67 


-1. 


.03 


-0.63 


47641. 


8367 


0. 


2486 


-116.31 


59.25 





.81 


-1.64 


47643. 


7878 





.7150 


60.26 


-113.94 


1. 


99 


1.53 


47644. 


8166 


0. 


9609 


-5.58 


-47.03 


1. 


20 


3.03 


47645. 


8002 





1960 


-111.58 


55.62 





.49 


-0.20 


47662. 


6849 


0. 


2320 


-117.86 


60.36 


-1. 


30 


0.03 


47664 


7592 


0. 


7279 


58.48 


-117.41 


-1. 


.07 


-0.66 


47674. 


8228 


0. 


.1334 


-93.47 


39.12 





.89 


1.12 


47675. 


.7352 


0. 


.3515 


-99.81 


43.29 


-0 


.14 


-0.06 


47689. 


7883 





.7107 


56.82 


-117.76 


-0 


.89 


-2.85 


47691. 


.7534 


0. 


1804 


-110.29 


53.39 


-1. 


.52 


0.89 


47691. 


9235 


0. 


.2211 


-116.21 


59.69 


-0 


.55 


0.26 


47693. 


8686 





6860 


54.15 


-110.29 





82 


0.21 


47693 


9378 





.7025 


56.23 


-112.40 


-0 


.26 


1.28 


47700. 


6733 


0. 


3126 


-109.29 


55.41 


1. 


.06 


1.32 


47723. 


6572 





8065 


56.53 


-111.37 


1. 


66 


0.68 


47725. 


6208 


0. 


2759 


-117.06 


61.70 


-1. 


.10 


1.98 


47727. 


6324 





7567 


58.66 


-118.05 


-1. 


.67 


-0.51 


47729. 


6085 


0. 


2291 


-116.36 


60.66 





00 


0.53 


47731. 


6271 


0. 


.7116 


57.97 


-113.55 





.14 


1.48 


47775. 


5008 


0. 


1989 


-112.54 


57.69 





.04 


1.36 


A 7777 


5 180 


0. 


.6811 


51.47 


-110.36 


-0 


.74 


-0.99 


47779. 


5228 


0. 


1603 


-103.69 


47.89 


-0 


30 


0.81 


47788. 


5003 





3062 


-111.36 


56.91 





29 


1.52 


48194. 


4589 


0. 


.3441 


-99.98 


43.14 


2 


.09 


-2.62 


48405. 


.7132 


0. 


8410 


46.85 


-104.96 





56 


-1.54 


48409. 


.7815 


0. 


.8135 


54.04 


-111.28 





60 


-0.67 


48429. 


7546 


0. 


5877 


17.51 


-73.70 


-0 


61 


1.40 


48432. 


7522 


0. 


3042 


-113.49 


54.28 


-1. 


.47 


-1.48 


48462. 


.5915 


0. 


.4369 


-63.20 


5.35 


-0 


.54 


-0.78 


48466. 


5983 


0. 


3946 


-81.02 


25.90 


1. 


.92 


-0.62 


48485. 


.5754 


0. 


9308 


8.77 


-65.98 


-0 


27 


-0.01 


48491. 


5446 


0. 


3576 


-98.40 


39.30 


-0 


.81 


-1.96 


48499. 


5903 


0. 


2808 


-114.65 


57.64 





.82 


-1.59 


48501. 


.5840 


0. 


.7574 


59.46 


-117.06 


-0 


.85 


0.46 



List of Objects 

'AD Boo' on pagej2 
'VZHya' on paged 
'WZ Oph' on pagejl 
'EWOri' on pageQ 
'UXMen' on page]! 
'HSHya' on page [3 
'V505Per' on pagey 
'ITCas' on page! 
'V1143Cyg' on paged 
'CDTau' on pagejl 
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'BWAqr' on pagei 



